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COMPUTER PROGRAMS FOR ESTIMATING AIRCRAFT TAKEOFF 
PERFORMANCE IN THREE-DIMENSIONAL SPACE 

Jeff Bowles 

Ames Research Center, NASA 
Moffett Field, Calif. 94035 

ABSTRACT 

A set of computer programs has been developed to estimate the takeoff 
and initial climb-out maneuver of a given aircraft in three-dimensional space. 
The program is applicable to conventional, vectored lift and power-lift con- 
cept aircraft. The aircraft is treated as a point mass flying over a flat 
earth with no side slip, and the rotational dynamics have been neglected. 

The required input is described and a sample case presented. 
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INTRODUCTION 


A set of computer programs has been developed to estimate the takeoff and 
initial climb-out maneuver of a given aircraft in three-dimensional space. 

The program is applicable to conventional, vectored lift and power-lift con- 
cept aircraft. The aircraft is treated as a point mass flying over a flat 
earth with zero sideslip, and the rotational dynamics have been neglected. 

The user is required to provide two subroutines which compute the total 
force coefficients along and normal to the flight path, and determine various 
required engine characteristics. 

This report describes the various subroutines and the required input, the 
equations used, and the computational techniques involved. Also included is 
a sample case and a listing of the program. 


NOTATION 


Symbol 

Fortran name 


a t 

T(6), S(10) 

acceleration along flight path (m/sec 2 , ft/sec 2 ) 

C L 

CL 

aircraft lift coefficient 

C D 

CD 

aircraft drag coefficient 

c x 

CX 

force coefficient along flight path 

C Y 

CY 

force coefficient normal to flight path and in 



plane of symmetry of the aircraft 

ENP 

ENP 

number of engines 

g 

G 

acceleration due to gravity (m/sec 2 , ft/sec 2 ) 

h 

HABS, S(7) 

altitude (m, ft) 

i 

w 

EYEW, EYEWNG 

incidence of wing (deg) 

LF 

XLF 

load factor 

q 


dynamic pressure (N/m 2 , lb/ft 2 ) 
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R/C 

S (15) , ROC, RTCL 

rate of climb (m/sec, m/min, ft/sec, ft/min) 

R/C . 
min 

ROCMIN 

minimum rate of climb during turning flight 



(m/min, ft/min) 

S 

SW, SWING 

reference wing area (m 2 , ft 2 ) 

T 

THRUST 

thrust, net or gross, per engine (N, lb) 

V 

T(4) , S(4) , VEL 

aircraft velocity (m/sec, ft/sec) 

V R 

VR 

rotation speed, EAS (m/sec, knots) 

W 

W, WG 

aircraft weight (N, lb) 

W f 

WF, WFUEL 

fuel flow (N/hr, lb/hr) 

X w 


wind axis coordinate, tangent to flight path (m. 

X 

S(8) 

earth fixed coordinate, along runway (m, ft) 

y w 


wind axis coordinate, perpendicular to x w and z w 



(m, ft) 

Y 

SC9) 

earth fixed coordinate, perpendicular to X and Z 



(m, ft) 

z w 


wind axis coordinate, perpendicular to x w and in 



plane of symmetry (m, ft) 

Z 

-S(7) 

earth fixed coordinate, normal to earth surface 



(m, ft) 

a 

ALPHA 

angle of attack (deg) 

Y 

S(5) , GAMMA 

flight path angle (rad, deg) 

6 f 

DELFD 

flap deflection (deg) 

6 s 

DELSPL 

spoiler deflection (deg) 

0 

THETAF 

pitch altitude (fuselage angle) (deg) 

y 

MU 

rolling coefficient of friction 

V 

ANGLE 

vectored thrust deflection angle (deg) 
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p RHO air density (kg/m 3 , slugs/ft 3 ) 

4, PHI roll angle, right wing down positive (deg) 

S(6) heading angle (deg) 

SUBROUTINE TAKOFF 

The subroutine TAKOFF simulates the takeoff and initial climb out maneuver 
of a given aircraft in three-dimensional space. The program is applicable to 
conventional, vectored lift and powered lift aircraft. The aircraft is treated 
as a point mass and the rotational dynamics of the aircraft are neglected. 

This simplification necessitates an estimation of all rotational rates involved. 
These rates are either input by the user or are approximated by a finite dif- 
ference form. In addition, the following assumptions are made: 

- flat earth 

- constant acceleration of gravity 

- zero sideslip angle 

The final assumption implies that the velocity vector and the resultant 
aerodynamic forces are contained in the plane of symmetry (ref. 1) . 

The takeoff maneuver is divided into four basic segments: ground roll 

and rotation, liftoff and initial segment climb, acceleration to final climb 
speed at a constant rate of climb, and finally, the pullup maneuver to 
establish the final climb speed. Provisions in the program are made for 
gear retraction, flap retraction, changing of vectored thrust angle and power 
setting, and changes in heading angle. 

The ground roll is made at a specified power setting and flap deflection. 
When the rotation speed is reached, the aircraft is "rotated" by increasing 
the angle of attack linearly with time until liftoff occurs or the tail scrape 
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angle is reached. If the latter occurs, the ground roll is continued with 
the fuselage angle equal to the tail scrape angle. 

The flight path control is obtained by monitoring four dynamic variables 
— acceleration along the flight path, load factor, fuselage angle (pitch 
attitude) and rate of climb. The aircraft is not permitted to decelerate and 
the load factor and fuselage angle are restricted to values less than or 
equal to a specified maximum value. If any of these conditions are violated, 
the angle of attack is reduced until all constraints are satisfied. During 
turning flight, if the rate of climb is less than a specified minimum value, 
the roll angle is reduced until the time rate of change of the rate of climb 
is non-negative. In addition, if the time rate of change of the flight path 
angle is less than -1.0 deg/sec, the roll angle is also reduced. 

Once a specified altitude is attained, called the maneuvering altitude, 
the aircraft is pitched down by a reduction in angle of attack until a specified 
rate of climb is obtained. The aircraft then accelerates at this rate of 
climb until the desired final climb speed is reached. 

When the final climb speed is attained, the pullup maneuver is executed 
in order to bring the aircraft to a zero rate of acceleration along the flight 
path. This maneuver is accomplished by increasing the angle of attack and 
pulling a load factor of 1.20, which will result in an increase in the rate of 
climb to a final value at the desired climb speed. It may be necessary to 
throttle the engines in order to maintain the desired constant climb speed 
subject to the fuselage angle restriction. 

Program Inputs 

The inputs to subroutine TAKOFF are through the argument list, input by 

NAMELIST and common blocks /UNIV/ and /AERO/. Either metric or English units 
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may be used in the program. On the first input data card, starting in col. 1, 
the word METRIC or ENGLISH should appear, depending on the user's choice of 
units. 

The call to TAKOFF is as follows: CALL TAKOFF (INPC, IDCN, WGROSS, 

SWING, XENG, VR, VEND) where 

INPC - program control = 1 - input data loaded 

- 2 - program executed 
= 3 - data input and program executed 
IDCN - print control = 1 - no print out 

= 9 - print out 

WGROSS - aircraft gross weight (N, lb) 

SWING - reference wing area (m 2 , ft 2 ) 

XENG - number of engines 

VR - rotation speed (m/sec, knots) 

VEND - final climb speed (m/sec, knots) 

All speeds are indicated in air speeds. 

There are three namelist inputs to TAKOFF, /NAM1/, /NAM2/ and NAM3/ . 

The namelist /NAM1/ input variables are as follows: 

CDGEAR - drag increment due to gear 
DFLPDT - flap retraction rate (deg/ sec) 

DTABS - temperature increment above standard temperature (°C, F) 

DTGR - time required to retract gear (sec) 

DTPDWN - throttling down rate (percent/sec) 

DTPUP - throttling advance rate (percent /sec) 

DTVECT - vectored thrust angle reduction rate (deg/sec) 

EYEWNG - wing incidence angle (deg) 
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HAPT - airport altitude (m, ft) 

HDT - obstacle height (m, ft) 

HGR - altitude at which gear retraction is started (n, ft) 

HMAN - maneuvering altitude (m, ft) 

HMAX - takeoff termination altitude (m, ft) 

UM - rolling coefficient of friction 
NPAGE - number of lines printed per page 
PMARG - pullup speed margin 

ROCMIN - minimum rate of climb during turning flight (m/min, fpm) 

ROLLMX - maximum allowable roll angle (deg) 

ROLRAT - roll rate (deg/ sec) 

RTCL - rate of climb during accelerate segment (m/min, fpm) 

THTFLY - maximum allowable fuselage angle while airborne (deg) 

THTSCP - tail scrape angle (deg) 

XLFMAX - maximum allowable load factor 

The user may input all, some, or none of the above input variables. The 
default values of these input variables are listed below: 

GDGEAR =0.0, DFLPDT =3.0 deg/sec, DTABS = 0.0° F, 

DTGR =5.0 sec, DTPDWN = 5.0 percent/sec, 

DTPUP = 6.0 percent/sec, DTVECT = 10 deg/sec, 

EYEWNG = 1.0 deg, HAPT = 0.0 ft, HDT = 35 ft, 

HGR =25.0 ft, HMAN = 1000 ft, UM = 0.02, 

PMARG = 0.04, ROCMDN = 250 fpm, ROLLMX =15.0 deg, 

ROLRAT = 5.0 deg/sec, RTCL = 750 fpm, 

THTFLY =15.0 deg, THTSCP = 10 deg, XLFMAX = 1.15. 

Note that all default values contained in the program are in English units. 
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If the default value of CDGEAR is used, the program will calculate, based 
on an empirical formula, a value for the gear drag as a function of gross 
weight and wing area. 

The second set of namelist variables, /NAM2/ , constitute the flap, 
throttle and vectored thrust schedules. These are tables that manage the 
flap setting, power setting and vectored thrust angle as a function of the 
aircraft speed and altitude. These variables are arrays of dimension 5. 

XDELFD (I) - flap deflection (deg) 

XHFLAP (I) - flap retraction altitude (m, ft) 

SVFLAP (I) - flap retraction speed (m/sec, knots) 

XPOWER (I) - power setting 

XHPWR (I) - power setting change altitude (m, ft) 

XVPWR (I) - power setting change speed (m/sec, knots) 

XNV (I) - vectored thrust angle (deg) 

XHVECT (I) - vectored thrust angle change altitude (m, ft) 

XWECT (I) - vectored thrust angle change speed (m/sec, knots) 

All altitudes are absolute altitudes and all speeds are indicated air 
speeds. 

These schedules are constructed as follows: If the speed or altitude of 

the aircraft is equal to, say, XVFLAP (I) or XHFLAP (I), respectively, then 
the flaps are retracted at the rate DFLPDT to the value XDELFD (I). The power 
setting and vectored thrust angle management work in a similar manner. The 
power setting may either be increased or decreased. The flap setting and vec- 
tored thrust angle setting can only be reduced with altitude or speed. The 
values of XDELFD (1), XNU (1) and XPOWER (1) are all for the ground run. The 
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user is permitted four changes in flap, power, vectored thrust angle settings 
during the airborne portion of the takeoff. 

The default values for /NAM2/ are as follows: 

- 100 percent throttle throughout takeoff 

- 0 degrees vectored thrust 

- 15.0 degrees flaps during ground roll, retracted to 5.0 degrees 
at 250-ft altitude, retraction to 2,0 degrees at 200 knots, and 
finally, complete retraction at 210 knots. 

Again, the user may choose to use all, some, or none of the above 
schedule values. Note that the default values contained in the program are in 
English units. No changes to any of these settings are allowed during the 
pullup maneuver. 

The final set of namelist variables, /NAM3/ , define the departure 
headings as functions of absolute altitude and ground distance from the start 
of takeoff roll point. The heading angle, with values - 180 < ip < 180, is 
positive for right turns proceeding along the flight path. These input 
variables are arrays of dimension 5. 

XHEAD (I) - flight heading (deg) 

XHHEAD (I) - heading change altitude (m, ft) 

XRANGE (I) - heading change ground distance (km, n.mi.) 

The departure heading schedule works in a similar fashion to the flap, 
power and vectored thrust angle setting schedules. If the absolute altitude 
or ground distance from the starting point of the takeoff roll is equal to 
XHHEAD (I) or XRANGE (I), respectively, the aircraft begins to turn to a 
heading value of XHEAD (I). The runway heading is defined to be a heading 
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angle of 0 degrees. Changes in aircraft heading are accomplished by increas- 
ing or decreasing (for right or left turns, respectively) the roll angle at a 
rate equal to ROLRAT. The absolute value of the roll angle is restricted to 
a maximum value of ROLLMX. The roll-out maneuver to establish the required 
heading is performed by rolling the aircraft back from the banked attitude to 
zero degrees roll (wings level) at a time such that when the wings are level, 
the aircraft is on the desired heading. The roll rate for the rollout maneuver 
is also equal to a value of ROLRAT. The default values for the heading 
schedule is for a straight out departure (no turns) . Four changes in heading 
angle are permitted during the takeoff. 

Program Output 

The program output consists of a time history of several aircraft and 
flight path parameters. The output will be in meters or English units, 
depending on the choice of the user. See the sample listing presented in 
Appendix B. The output variables are as follows: 

TIME - time from start of takeoff roll (sec) 

X DIST - ground track distance along the earth fixed X coordinate 

(m, ft) 

Y DIST - ground track distance along the earth fixed Y coordinate 

(m, ft) 

ALT - aircraft altitude (m, ft) 

TAS - true airspeed along flight path (m/sec, knots) 

EAS - indicated airspeed (m/sec, knots) 

MACH NO - Mach number 

ACCEL - acceleration along flight path (m/sec 2 , ft/sec 2 ) 

CL - aerodynamic lift coefficient 
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CD - aerodynamic drag coefficient 

ALPHA - angle of attack (deg) 

GAMMA - flight path angle (deg) 

R/C - rate of climb (m/min, fpm) 

LOAD FACTOR - load factor 

THRUST - total thrust, net or gross (N, lb) 

FUS. ANG. - fuselage pitch angle (deg) 

ROLL ANGLE - roll angle (deg) 

HEADNG - heading angle (deg) 

In addition, the user may also obtain the following values through the 
common block /EXCHNG/ : 

SROLL - distance to liftoff (m, ft) 

S35 - track distance to obstacle height (m, ft) 

V35 - speed (EAS) at obstacle height (m/sec, knots) 

The program will terminate normally when the end speed is reached 
(VEND) or when the maximum specified altitude (HMAX) is attained. Abnormal 
termination will occur under several conditions: 

- flight path constraints cannot be met by further reduction 

in angle of attack 

- aircraft cannot accelerate at input rate of climb (RTCL) 

- aircraft altitude becomes negative 

- ground track distance in ±X or ±Y direction is greater than 

10 n.mi. 

- ground run exceeds 90 sec 

- elapsed time greater than 300 sec 
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For further definitions and explanations refer to the listing of TAKOFF 
and supporting subroutines contained in Appendix C, and the sample case pre- 
sented in Appendix B. 

Subroutines ARODYN and ENGINE 

The takeoff subroutine described above requires the user to provide two 
subroutines to compute total force coefficients and determine various required 
engine characteristics (e.g., thrust and fuel flow per engine). The format 
and structure of these subroutines is left to the discretion of the user. The 
units used in these subroutines should be the same as those of the input data. 


I . Subroutine ARODYN 

This subroutine computes the total force coefficients along the flight 
path and normal to the flight path in the plane of symmetry as a function of 
angle of attack and thrust. A force coefficient in a particular direction 
e s is defined to be the sum of all aerodynamic and propulsion system forces 
in that particular direction, divided by the dynamic pressure times the wing 
area. 

A e e ■ EF 

C k 

s qS 

The transfer of the various computer variable values to and from sub- 
routine ARODYN is through labeled common blocks /UNIV/ and /AERO/. Of primary 
concern is the common block /AERO/ : 

COMMON /AERO/ VEL, QS, HABS, THRUST, TVECT, ANGLE, DELFD, DELSPL, ALPHA, 
CS, CY, CL, CD, RHO, GRCD, I FAST 

The input variables from TAKOFF are: 

VEL - aircraft velocity along flight path (m/sec, ft/sec) 
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QS 

- dynamic pressure times wing area 

(N, 

lbs) 

HABS 

- absolute altitude of aircraft (m. 

ft) 


THRUST 

- thrust (net or gross) per engine 

(N, 

lbs) 

TVECT 

- total vectored thrust (N, lbs) 



ANGLE 

- angle of vectored thrust relative 

to 

aircraft center line. 


positive down (deg) 



DELFD 

- flap deflection (deg) 



DELSPL 

- spoiler deflection (deg) 



ALPHA 

- angle of attack (deg) 



RHO 

- air density (kg/m 3 , slugs/ft 3 ) 



GRCD 

- drag increment due to gear 




The return from ARODYN should be: 

CX - total force coefficient along flight path 

CY - total force coefficient normal to flight path in plane of 

symmetry 

The output variables CL and CD are provided to the user as a means to 
distinguish between pure aerodynamic coefficients and total force coefficients. 
The output variables CL and CD are printed out in the time history, but are 
not used in the actual calculations. If desired, in subroutine ARODYN, CL and 
CD may be directly equated to CY and CX, respectively. 

There is a certain amount of redundancy among some of the input variables. 
The user may utilize only those variables he desires and disregard the others. 
Due to the wide range of velocities encountered during the takeoff, there will 
be a correspondingly large variation in the magnitude of the force coefficients 
which must be accommodated in subroutine ARODYN. 
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II. Subroutine ENGINE 


This subroutine provides the various propulsion data to subroutines 
TAKOFF and ARODYN. 

The inputs to subroutine ENGINE are through the argument list and labeled 
common blocks /AERO/ and /UNIV/ . 

The call to ENGINE is as follows: 

CALL ENGINE (ALT, DTABS, EN, PWRSET, WFUEL, KENG) 

where 

ALT - aircraft altitude (m, ft) 

DTABS - temperature increment above standard temperature (°C, °F) 

EN - aircraft Mach number 

PWRSET - power setting (see below) 

WFUEL - fuel flow (N/hr, lbs/hr) 

KENG - engine control parameter = 0 


The variable PWRSET is defined to be: 


PWRSET = 


net thrust 

net thrust available 


and is the parameter used in controlling the thrust level. It is used for 
power setting management during the takeoff. 

The user may choose to work with either the gross thrust per engine or the 
net thrust per engine, provided he uses the variable THRUST properly in the 
calculation of the total force coefficients. For example, when using gross 
thrust per engine, the ram drag must be included in the total summation of 
forces. If the gross thrust vector and ram drag vector can be considered 
col linear, the user may choose instead to work simply with the net thrust. 

Refer to the sample case presented for an illustration of subroutines 
ARODYN and ENGINE. 
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APPENDIX A 


EQUATIONS 

1. Equation of motion during ground roll 

dV/dt = (g/W)[-Wp + qS(C y y - ] 

2. Equation of motion along flight path 

dV/dt = (g/W) (-C x qS - W sin y) 

3. Equation of motion normal to flight path and in the plane of symmetry 

d'p/dt sin <J> cos y + dy/dt cos = (g/WV) (C y qS - W cos y cos <j>) 

4. Equation of motion normal to the flight path and normal to plane of 
symmetry 

-di|»/dt cos <f> cos y + dy/dt sin 4 = (g/WV) sin <t> cos y 

where 

g = gravity constant 
W = aircraft weight 
q = dynamic pressure 
S = wing area 
y = flight path angle 
iJj = heading angle 
<f> = roll angle 

V = aircraft velocity 

C x = total force coefficient along flight path 

C y = total force coefficient normal to flight path and in plane of 
symmetry 
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The coordinate systems used are presented in figure 1. The XYZ is the 

right handed earth fixed coordinate system, with the X-axis aligned with the 

runway. The Z axis is vertical and positive downward. The wind axis system 

is defined as follows: the x w axis is tangent to the flight path and positive 

in the direction of flight; the z is normal to the x axis, in the plane of 

w w 

symmetry of the aircraft, and is positive downward in level flight; the y w 

axis is normal to the x and z axis in the right handed sense. The 

w w 

x, y, z coordinate system is the translation of the XYZ axis system to the 
location of the point mass representation of the aircraft (ref. 2) . 

In order to make the system of equations of motion more amenable to 
numerical integration, the equations are manipulated in order to obtain 
explicit relations for the time rates of change of the velocity, flight path 
angle and heading angle. 

The equation for dV/dt is already in the desired form. Note that the 
acceleration along the flight path is independent of the roll attitude. To 
obtain an expression for dy/dt alone, equation 3 is multiplied by cos <j>, 
equation 4 multiplied by sin <f> and the resulting equations subtracted to 
give: 

5. dy/dt = [g/(WV)](C y qS cos $ - W cos y) 

An expression for di^/dt alone is obtained in a similar manner: 

6. di(i/dt = [g/(WV cos y)] C y qS sin 4 

The system of equations 2, 5, and 6 are then numerically integrated 
using the Adams-Moulton fixed step-size method. 

7 . Load factor 
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8. Constant rate of climb equation 

Rate of climb = ROC ~ V sin y 

For ROC to be constant with time, 

dROC _ 
dt = U 

or, 

dROC d . dV . „ dy 

It" = dt (V Sin y) * dt sin Y + V cos Y dt 0 

Substituting for terms dV/dt and dy/dt from equations 2 and 5, and 

simplifying: 

qS(C cos y cos 4> - C sin y) - W = 0 

y * 

9, Rotational rate approximations by finite difference 

6 = y + a - i 
' w 

where 

6 - pitch attitude (fuselage angle) 

y - flight path angle 
a - angle of attack 
i - incidence of wing 

Differentiating with respect to time we obtain: 

d8 _ da dy 
dt ~ dt + dt 

(jv 

where jj: is given by equation 5. 

da . 

^ is approximated by the finite difference form: 

dt - ^°now a past^ At 
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where 


a = current value of the angle of attack 
now 

°past = P rev ^ ous vaque angle of attack 
it = integration time interval 
10. Roll-out control equation 

As the desired heading angle ip£ is approached with the aircraft 
banked at some angle of roll the roll angle is reduced at the rate of roll 
ROLRAT to zero in such a way that when the wings are level (implying 
diji/dt = 0), iJj = ip f . To perform this roll-out maneuver, an open loop type 
control procedure is used. The problem is to determine at what heading 
angle ip the roll-out should be initiated. 

From equation 6 

d* sqsc gqsc 

dt - WV cos y S '*' n ^ WV cos y ^ 
for moderate angles of bank. 

Using a finite difference form approximation for di^/dt and the defini- 
tions of figure 2, 


d^ 


»f ~ * _ A* m gqSC y 


t^ WV cos y ave 


dt At 

It is desired that the time average of <f> over the time interval 
0 < t’ < equals where 

I clV C 

= 4> t + t T » ^ being constant 

therefore 

-t. 


<j> t. 
Y ave f 


■ /v **■)«• 
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or, upon integration, 


10a. 


Ajp wv cos Y t + 1_ d<j> 2 

t £ gqSC y r f " Vf 2 dt T f 


Now, for ^(t^) = 0, 



Solving for t £ and substituting in equation 10a 


A* = 


<t> T 2 gqsc y 

2(d<J)/dt)WV cos y 


The value of A^ is monitored during turning flight, and whenever 
l+f ~ S | A^ | , the roll-out maneuver is begun. This estimate of Aip is 

not exact, since the velocity V and the flight path angle y will change 


over the time period t £ , but for moderate roll angles and roll rates, t £ 
will be small, and hence changes in V and y correspondingly small. 
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time t f 

t'=0 - start of rollout maneuver 
4(0) = » ♦ Uf) = 0 

Figure 2.- Rollout maneuver parameters 


t l 
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APPENDIX B 


SAMPLE CASE 

Shown below is an example of the input, calling format, subroutines 
ARODYN and ENGINE, and the print out obtained from the takeoff program. 

The main calling program TEST1 is set up to do the takeoff of a 
Boeing 727-200. The required common blocks are shown, but others may be added 
if needed. This example was done in English units. 

The input was as follows: 

ENGLISH 

$NAM1 NPAGE = 48, RTCL = 550., 

THTFLY = 20., HMAN = 2000., 

ROLLMX =30., R0CMIN - 500. $END 
$NAM2 XPOWER(2) = 0.75, XHPWR(2) = 750., 

XP0WER(3) = 0.95, XHPWR(3) = 1750. $END 
$NAM3 XHEAD(l) = 45., XHHEAD(l) = 800., 

XHEAD(2) = -15., XHHEAD (2) = 2250. $END 

Subroutine ARODYN calculates the lift and drag coefficients of the 727-200 
as a function of angle of attack, flap and spoiler deflection. The increments 
of lift and drag due to flaps is determined by a table look-up format. Once 
the lift and drag coefficients are computed, the thrust components, normalized 
by dynamic pressure times wing area (QS) , are added in to determine the total 
force coefficients CX and CY. 

Subroutine ENGINE computes the thrust and fuel flow of the JT8D engine, 
based on a simplified model. The thrust lapse is assumed to be linear with 
Mach number, and the fuel flow assumed linear with power setting. 
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This particular run was made on the Lawrence Berkeley Laboratory 
CDC 7600, requiring a field length of 41700 words to load and 2.47 sec to 
execute. 
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TEST t 


3 

5 

b 

7 

10 

12 

14 

22 


PROGRAM TEST 1 < INPUT, OUTPUT, TAPES* INPUT, TAPE6«0UTPUT) 

COMMON /UN I V / NPC # NSC ,I0C #H , ST ,R ,W , 

IMF ,fM ,VMO , F. H M 0 , ALPHL 0 , C L AUPH, S M ,AR ,8 # 

2tYtw ,FNP , T A ,*G ,wGS , K WR I TE , DL HC4 

3, KS 1 7E 

common /aero/ vfl,gs, nabs, thrust, tvect, angle, delfd»dclspl, alpha, 
9CX#CY,CL,r0,RH0iGRCr>, IF AST 
wG = 172000* 

S*ING c 1720 . 
fNP c J.O 
OELSPL = 0.0 
RHO » 0.0023 
S* « SWING 
H ■ WG 

CALL TAKOFF (3,9, WG,SW, 3.0, 135,0, 250, 1 
END 


rsj 

4 ^ 


PROGRAM LENGTH INCLUDING 1/0 BUFFERS 

01116 

FUNCTION ASSlGNMFfcTS 
STATEMENT ASSIGNMENTS 
BLOCK NAMES ANO LENGTHS 


UNIV • 000030/01 

Af RO 

• 000020/02 





VARIABLE assignmfnts 




sw • 000015/01 swing • 000047 


OELSPL - 000007/0? 
MG - 000023/01 

Enp 

- 000021/01 

RHO * 

000OJS/02 

w 


START OF CONSTANTS- 

0OOO2S 

i 

1 

<r, 

a. 

UJ 

OOOOqS 

INPIRECTS- 

000047 


7600 COMPILATION — 

RUN76 

LEVEL RB 

74/07/15, 





ROUTINE COMPILES IN 04<JQi>0 


• 000006/01 



ARODYN 


N> 

Ln 


SUBROUTINE AR0DV\ 

REAL Mu 

CGHMO* /UNI V/ nPC ,N$C iIOC ,M ,9T ,R ,W ' 

IMP ,Fm ,VhO ,E*M0 , ALP h lO#Cl*LPH# SW »A« rB f 

2EYEW ,F.NP , T A *WG #WGS . KrfRITE* DLMC4 

«g£ /AfRO/ VEL,0S»HA6S, THRUST, TVPCT, ANGLE # DELFD » DFLSPL » ALPHA* 
9CX*CY,CL»n>#RH0*GRCP,IFAST 
DIMENSION XDElFri(6),XOfcLCLtA)»XOELCOC6).XVDELG(6) 

DAI A XDE LF 0/ 0 * 0 r S . 0 , 10, ,15, *20, ,25,/ 

DATA XDELCL/Q, 0,0,1 8b, 0,347, 0,482,0, 60, 0,702/ 

data xdelcc/o ,q» , oiae,,0295,,o45i,,0607, t o837/ 

DATA XVOEI 6/ 1,0 ,,995,, 990, ,980, ,970,,955/ 

CALL lTRLN(XOFLFO,XDElCD*DELFDf PEl.C0F,6) 

5 CALL lT*LMXDFLFD,XDELCL,PEt FD,DELCLF,6I 

11 CALL I TRl^CXPELFD. XVDEL6,0EIFD, VDFL6f 6) 

17 S A5 5 o ^ Q I ei 

20 CLAIPH s 4,5 

22 ALPHLO = -1*5 

23 SA7 * 0,0548 

25 SIGMA * 0,8 

26 DCLSPl = «,3lAr(0tl SPL/90,) 

30 DCP5PL = O,!2a(DEI$PL/90 J 

32 CL * CLALPH* c ALPHA - AlPHLO)* ,017453 ♦ DELCIF 

^0 CD ' SA5 f D r>ELCOF t ( SA 7 / VDEL6 ) * (CL * $ I GM A*DLLC LF ) * *2 + GRCD 

47 CO A CD + ntOSPL 

51 1 AlPHX 3 ALPHA * ,0174533 

53 IF ( 05 ,E0, 0.0>GS = 0,1 

55 CX = FO - THRU5T*FNP*C0S(ALPHX) /QS 

63 CY * Cl. * THRtJST*FNP*SlN(AlPHX)/Q8 

72 RETURN 

72 END 


subprogram length 


00167 


FUNCTION ASSIGNMENTS 

statement ASSIGNMENTS 
1 - 00005? 

Block names and lengths 

UNIV - 000030/01 aero - 000020/02 


VARIABLE ASSIGNMENTS 

ALPHA . 000010/02 ALPHLO - 00 OC 1 3/0 I AlPHX 

CX - 000011/02 CY - 000012/0? DCHSPL 

OELFO - 000006/0? DLL SPl. - 000007/02 ENP 

SA5 - 000161 SA7 - 000162 SIGmA 

XDELCL - 000134 XDELFD - 00012b XVDFL6 


000166 CD 

000165 DCLSPL 

000021/01 GRCD 
000163 THRUST 

OOOlSO 


START OF CnNST ANTS-0000 75 TEMPS--00 0 1 1 2 I wD I REC TS" 0 00 1 25 


000014/02 CL • 000013/02 CLALPH 

000164 DEL CDF • 0O01S6 DfcLCLF 

000016/02 Nil • 000125 OS 

000003/02 V0EL6 - 000160 XDElCD 


000014/01 

000157 

000001/02 

000102 



ENGINE 


SUBROUTINE Engine r At T # DTAHSiEN, PNRSF.T#*FuELrKENG> 

COMMON /AFRO/ VEL,0S,HA8S f THRUST, TVECT, ANGLE, Of LF P, OElSPL# ALPHA, 
9CX,CY,CL,CD,RN0,GRCO, If AST 


11 


IFCKENG # F0, 1 ) GO TO 10 
TO * 14000. 

12 


THRUST * f T f! - 6,O*ENM100 # )*PWRSET 

16 

11 

WFUEL * THRUST*0,63*PWftSET 

20 


RETURN 

21 

10 

PWRSET * THRtlST/fTO - 6,0* tM *1100,) 

25 


GO TO 11 

26 


END 


SUBPROGRAM length 
00043 

FUNCTION ASSIGNMENTS 

STATEMFNT ASSIGNMENTS 

10 - 000022 11 - 000017 

BLOCK names and LENGTHS 

AERO - OOC 020/0 1 

VARIABLE ASSIGNMENTS 

THRUST - 000003/01 TO - 0 0Q042 

START OF CONSTANTS-00003 I TtwPS — OOOO 36 I NO I RFC IS- 0 00 0 4 2 

7600 COMPILATION — PUN76 lEVEl 0* 74/07/15, 

ROUTINE COMPILES IN 044000 



»• INUjt* to T*« OFf • ALTITUDE • 


0,0 tEHPERATU*E ■ 


54,0 DEC 


A/C CHARACTERISTICS 

GROSS RAMI* MT, ■ 172000 MING AREA ■ 1720 

MING LOADING • 100,0 THRU8T/NEIGHT ■ ,240 


STATIC SEA LEVEL THRUST • 14000 


A/C PARAMETERS# 

NO, ENGINES ■ 3,0 


CDOEAR • ,0297 


EYEMNG « 1,0 


TAIL SCRAPE ANGLE ■ 10,0 


plight path control parameters 

MAX LOAD PACTOR ■ 1,10 GEAR RETRACTION ALT, ■ 25,0 

MANEUVER ALT, » 2000 ACCELERATE RATE OF CLIMB • 


MAX PLOOR ANGLE ■ 20,0 
550 


parameter variation rates 

DAOT ■ 1,0 DPLPOT » 3,0 DTGR • 5,0 

DTOUP • 6,0 OTVECT * 10,0 


DTPOMN « 5,0 


POMER, VECTORED THRUST, AND PLAP SCHEDULES 


to 

"J 


THROTTLE/POMER SETTING 
PNR8ET 1.00 

• 75 


1,00 

1,00 

SPEED 

0,0 

0,0 

0,0 

999,0 

999,0 

ALTITUDE 

0 

750 

1T50 

0 

0 

VECTORED THRUST 

angle 

angle 

0,0 

0.0 

0,0 

0,0 

0,0 

SPEED 

0,0 

900,0 

900,0 

999,0 

999,0 

ALTITUDE 

0 

0 

0 

0 

0 

PLAP DIPLCCTION 

deled 

ANGLP 

15,0 

5,0 

2,0 

0,0 

0,0 

SPEED 

0,0 

0,0 

200,0 

210,0 

0,0 

ALTITUDE 

0 

250 

0 

0 

0 


ALL 8PEC0S ARE INDICATED AIR SPEEDS AND ALL ALTITUDES ARE ABSOLUTE ALTITUDES 


DEPARTURE HEADINGS 


RANGE 

100,0 

100,0 

100,0 

100,0 

100,0 

ALTITUDE 

800,0 

2250,0 

99999,0 

99999,0 

99999,0 

HEADING 

A5,0 

•15,0 

0,0 

0,0 

0,0 



TAKEOFF (ELEVATION ■ 


0 FT) 


TIME 

X DI 5 T 

V DIST 

ALT . 

TAS 

EAS 

MACH 

ACCEL 

CL 

CD 

ALPHA 

GANNA 

R/C 

LOAD 

THRUST 

FU8, 

ROLL 

HEA 0 N 6 

( SEC ) 

trtET ) 

( FEE T ) 

( FEET ) 

( KTS ) 

( KT 8 S 

NO . 

( FP «2) 



( DEO ) 

( DEC ) 

( FPM ) 

FACT 

( LBS ) 

ANC ; 

ANCLE 

( DEG ) 

0.0 

0.0 

0,0 

0,0 

0,0 

0.0 

0,000 

7,22 

,470 

,0902 

1,00 

0,00 

0,0 

0,00 

42000 

o.o 

0,0 

0,0 

>.# 

3.6 

o.o 

0.0 


«.3 

,006 

7,20 

,678 

,0902 

1,00 

0,00 

0,0 

0,00 

41605 

0,0 

0,0 

0,0 

2,0 

14,4 

0,0 

0,0 

6,5 

8,5 

,012 

T , 1 7 

,678 

,0902 

1.00 

0,00 

0.0 

0,00 

41750 

0,0 

0,0 

0,0 

1,0 

32.4 

0.0 

0,0 

12,# 

12,* 

,019 

M 4 

,670 

,0902 

1,00 

0,00 

0,0 

0,00 

41431 

o.o 

0,0 

0,0 

0,0 

57.5 

0,0 

0.0 

17,0 

>7.0 

,025 

7,10 

,678 

,0902 

1,00 

0,00 

0,0 

0,00 

41504 

0.0 

0,0 

0.0 

5,0 

69,7 

0,0 

0,0 

21,2 

21,2 

,031 

7,06 

,670 

,0962 

1,00 

0,00 

0,0 

0,00 

41379 

0.0 

0,0 

0,0 

*.0 

129,0 

0,0 

0,0 

25,3 

25.3 

,010 

7,02 

,678 

,0982 

1,00 

0,00 

0,0 

0,00 

41254 

0,0 

O.o 

0,0 

T.O 

175.3 

0,0 

0,0 

20,5 

20,5 

,044 

6,96 

,670 

,0982 

1,00 

0,00 

0,0 

0,00 

41130 

0,0 

0,0 

0,0 

8,0 

220.6 

0.0 

0,0 

11,6 

31,6 

,050 

6,?3 

,670 

,0902 

1,00 

0,00 

0,0 

0,00 

41006 

0,0 

0,0 

0,0 

0,0 

206.0 

0,0 

0,0 

17.7 

37.7 

,056 

6,66 

,670 

,0902 

1,00 

0,00 

0,0 

0,00 

40884 

0,0 

0.0 

0,0 

10,0 

355.9 

0,0 

0,0 

41 , S 

41.8 

,063 

6,63 

,670 

,0902 

1,00 

0,00 

0,0 

0,00 

40762 

0,0 

0,0 

0,0 

11,0 

429.0 

0,0 

0,0 

45.8 

45,6 

,069 

6,76 

,678 

,0902 

1,00 

0,00 

0,0 

0,00 

40642 

0.0 

0,0 

0,0 

12,0 

510.5 

0,0 

0,0 

40,8 

40,6 

.075 

6,72 

,670 

,0902 

1,00 

0,00 

0,0 

0,00 

4 Q 522 

0,0 

0.0 

0,0 

11,0 

597,9 

0.0 

0,0 

51,7 

55.7 

,061 

6,66 

,678 

,0902 

1,00 

0,00 

0,0 

0,00 

40403 

0,0 

0,0 

0,0 

H . 0 

692.0 

0.0 

0,0 

57,7 

57,7 

,067 

6,60 

,670 

,0902 

1,00 

0,00 

0,0 

0.00 

40266 

0,0 

0.0 

0,0 

15.0 

792,6 

0,0 

0,0 

61,5 

61,5 

,092 

6,53 

,670 

,0902 

1,00 

0,00 

0.0 

0.00 

40169 

0,0 

0,0 

0,0 

10,0 

099.0 

0.0 

0,0 

65,4 

65,4 

,090 

6,47 

,670 

,0962 

1,00 

0,00 

0,0 

0,00 

40054 

0,0 

0,0 

0,0 

l *,0 

10 t 3,5 

0,0 

0,0 

60,2 

60,2 

,104 

6,40 

,676 

,0982 

1,00 

0,00 

0,0 

0,00 

39940 

0,0 

0,0 

0,0 

1«.0 

1133,6 

0.0 

0,0 

73.0 

73.0 

,uo 

6,33 

.670 

,0902 

1.00 

0,00 

0,0 

0,00 

39027 

0,0 

0.0 

0,0 

1*,0 

1260,0 

0.0 

0,0 

76,7 

76.7 

,115 

6,25 

,678 

,0902 

1.00 

0,00 

0,0 

0,00 

39715 

0,0 

0.0 

0,0 

20,0 

1392.6 

0,0 

0,0 

80,4 

60.4 

,121 

6,10 

,670 

,0902 

1,00 

0,00 

0,0 

0,00 

39605 

0.0 

0.0 

0,0 

21,0 

1531.4 

0,0 

0,0 

84,0 

84,0 

,126 

6,10 

,670 

,0902 

1.00 

0,00 

0,0 

0.00 

39496 

o.o 

0,0 

0,0 

22,0 

1676,4 

0,0 

0,0 

87,6 

87.6 

.132 

6,03 

,676 

,0962 

1.00 

0,00 

0,0 

0,00 

39380 

o.o 

0.0 

0,0 

21,0 

1027.3 

0,0 

0,0 

01,1 

01,1 

,117 

5,95 

,676 

,0902 

1.00 

0,00 

0,0 

0,00 

39202 

0,0 

0,0 

0,0 

24,0 

1904,2 

0,0 

0,0 

Oil , 6 

04.6 

#143 

5,07 

,670 

,0902 

1,00 

0,00 

0,0 

0,00 

39177 

0,0 

0.0 

0,0 

25,0 

2147.0 

0,0 

0,0 

08,1 

08,1 

,146 

5,79 

,676 

,0902 

1,00 

0,00 

0,0 

0,00 

39073 

0,0 

0,0 

0.0 

20,0 

2315,6 

0,0 

0,0 

101,5 

101,5 

• 153 

5,71 

,670 

,0902 

1,00 

0,00 

0,0 

0,00 

30971 

0,0 

0,0 

0,0 

21,0 

2409,0 

0,0 

0,0 

104,0 

104,0 

,156 

5,62 

,676 

,0902 

1,00 

0,00 

o.o 

0,00 

30071 

0,0 

0,0 

0,0 

20,0 

2669,7 

0,0 

0,0 

108,2 

108,2 

.103 

5,54 

,678 

,0902 

1,00 

0,00 

0,0 

0,00 

30772 

0,0 

o.o 

0,0 

20,0 

2055.2 

0,0 

o.o 

111,4 

111,4 

,160 

5,46 

,676 

,0902 

1,00 

0,00 

0,0 

0,00 

30674 

0,0 

0,0 

0,0 

10,0 

3046,1 

0,0 

0,0 

114,6 

114,6 

,173 

5,37 

,670 

,0902 

1,00 

0,00 

0,0 

0,00 

38570 

0.0 

0.0 

0,0 

11,0 

3242,3 

0,0 

0,0 

117,8 

117,8 

,170 

5,29 

,678 

,0902 

1,00 

0,00 

0,0 

0,00 

36463 

0.0 

0,0 

0,0 

12,0 

3443,9 

o.o 

0,0 

120,0 

120,9 


5,20 

,670 

,0982 

1,00 

0,00 

0,0 

0,00 

30390 

0,0 

0.0 

0,0 

31.0 

3650.7 

0.0 

0,0 

123,0 

123.9 


5,12 

,670 

, 0* S 2 

1,00 

0,00 

0,0 

0,00 

30299 

0,0 

0,0 

0,0 

34.0 

3662,5 

0,0 

0,0 

126,0 

126,0 

.101 

5,03 

,678 

,0962 

1,00 

0,00 

0,0 

0,00 

30209 

o.o 

0.0 

0,0 

55,0 

4079,4 

0,0 

0,0 

120,0 

120,0 

,196 

4,94 

,670 

,0902 

1.00 

0,00 

0,0 

0.00 

36120 

0,0 

o.o 

0,0 

10,0 

4301.3 

0.0 

0,0 

112.8 

132.8 

,200 

4,06 

,670 

,0902 

1,00 

0,00 

0,0 

0,00 

38033 

0,0 

0.0 

0,0 

ROTATION (TIME • 

36,6 

AND TAfi 

1 * 135.1 

EAS « 

135,1) 











37,0 

4520.0 

0,0 

o.o 

135,6 

135,6 

,205 

4,76 

,694 

,0909 

1,20 

0,00 

0,0 

0,00 

37940 

.2 

0,0 

0,0 

30.0 

4754,4 

0,0 

0.0 

130,4 

130,4 

,209 

4,63 

,773 

,1020 

2,20 

0,00 

0,0 

0,00 

37064 

1.2 

0,0 

O.o 

34,0 

4045,5 

0,0 

0,0 

141,1 

141,1 

.313 

4,47 

.8*1 

.1074 

3,20 

0,00 

0,0 

0,00 

37703 

2,2 

o.o 

0,0 

40,0 

5210,1 

0,0 

0,0 

143,7 

143,7 

,217 

4,29 

,930 

,1127 

4,20 

0,00 

0,0 

0,00 

37705 

5,2 

0,0 

0,0 

41,0 

5400,0 

0.0 

o.o 

146,2 

146,2 

,221 

4,10 

1,006 

,1106 

5,20 

0,00 

0,0 

0,00 

37631 

0.2 

0,0 

0,0 

42,0 

5720,0 

0,0 

0,0 

140,6 

140,6 

,224 

3,00 

1.067 

,1252 

6,20 

0,00 

0,0 

0,00 

37560 

5.2 

0,0 

0.0 

41,0 

5002,6 

0,0 

0,0 

150.0 

150,6 

,220 

3,63 

1,165 

,1326 

7.20 

0,00 

0,0 

0,00 

37493 

6,2 

0,0 

0,0 

LIFTOFF 

(TIME « 

43,9 

DIST « 

6213,3 TAS 

• 152.7 EAS 

* 152 

:#n 









44,0 

6230,1 

0,0 

.0 

152,0 

152,9 

,231 

5,35 

1.244 

,1406 

6,20 

.01 

3.0 

1.02 

37430 

7,2 

0,0 

0,0 

45,0 

6400,0 

0.0 

.0 

154, T 

154,7 

,234 

2.70 

1,311 

,1479 

9,05 

.48 

126,2 

1,10 

37374 

6,3 

0,0 

0,0 

44,0 

6761,4 

0,0 

4.6 

156,2 

156,2 

1 236 

2*30 

1,287 

,1453 

0,75 

1.1 S 

310,6 

1,10 

37326 

0,9 

0.0 

0,0 

4 T ,0 

7026,2 

0,0 

11.5 

157,5 

157,5 

,238 

2,00 

1,267 

,1431 

8,50 

1.64 

511,7 

1.10 

37260 

6,3 

0,0 

0,0 

40,0 

7203,0 

0,0 

21.6 

158,6 

150,5 

,240 

1,62 

1,240 

.1410 

0,25 

2,52 

707,1 

1.10 

37255 

6 ,« 

0,0 

0,0 



Takeoff continued 


TIN! 

X OUT 

y out 

ALT. 

IFEEf) 

TA8 

EAS 

NACf* 

ACCEL 

CL 

(ICC) 

(fCCTj 

(FEET) 

(XT*) 

(XT*) 

NO, 

(pp*a> 


SCAR RETRACTION 

started 

AT 45.1 

•Ce, COMPLETE 

AT 

9 i.i tec 


DISTANCE TO OMTACLE HEIGHT ■ 7560. A SCREEN SFEED(EAS) « 159 

*4 

44,0 

7941,1 

0,0 

19,0 

194,9 

156,4 

,201 

1,** 

1,214 

96,0 

7410,7 

0,0 

91,7 

140,2 

160,1 

,*•* 

1,04 

1,228 

St ,0 

•101,0 

0,0 

71,* 

1*0, T 

160,6 

,*01 

,*T 

1,220 

91,0 

•171,4 

0,0 

49,1 

141,* 

161,0 

.*00 

,4« 

1,212 

91,0 

•441,1 

0,0 

121.0 

141,9 

161,2 

,*00 

.«* 

1,20* 

94,0 

•*14,4 

0.0 

191.* 

1*1,7 

161,3 

,*09 

. 0 * 

1,20* 

99,0 

4189,4 

0,0 

1*4,4 

1*1,7 

161,2 

,*09 

,00 

t,l«* 

94,0 

RL*R» 

4494,1 

RCTRACTCO 

0.0 
70 9,0 

21*;* 
DEG, IN 

141,7 
1,1 SEC 

161,1 

rt 

.209 

.00 

1,110 

97,0 

6717, 1 

0,0 

294,4 

1*1,7 

1*1,1 

,109 

.02 

1,102 

99,0 

6447,4 

O.o 

««*;* 

1*1,7 

1*1.1 

.*09 

.1* 

1,100 

96,0 

10244,4 

0,0 

129,9 

141,* 

1*1.1 

,*09 

,20 

1,041 

40,0 

10940,0 

O.o 

140,8 

1**,0 

1*1.* 

,I« 

.17 

1,072 

41,0 

10«t 1,9 

0,0 

1*4.1 

14*,* 

141.1 

.*«* 

.19 

1,12* 

6 t ,0 

11081.0 

0,0 

411.0 

14*,* 

141.2 

.*«* 

,00 

1,117 

49,0 

11194,4 

O.o 

470,4 

1**,* 

1*1,1 

,*0* 

.01 

1.041 

44,0 

11429,7 

0,0 

904,1 

14*,* 

141.0 

,*«4 

,01 

1,0*9 

49,0 

11*67,0 

0,0 

so7;s 

1**.* 

1*0,4 

,*04 

,00 

1,0*9 

44,0 

12140,1 

0,0 

9*9,4 

14*,* 

140,4 

,*04 

,00 

1,0*9 

47,0 

12414,4 

0,0 

424,4 

14*. * 

140,* 

,104 

.00 

1,0*5 

40,0 

1*710,4 

0,0 

442,8 

14*,* 

1*0,7 

.104 

.01 

1,0*9 

44,0 

1*442,2 

0,0 

701,1 

14*,* 

140,4 

,204 

.01 

1,0*5 

70,0 

rctord 

11*91,9 

0,0 

714,7 

14*. 1 

1*0.9 

. ,204 
,0 IE< 

,01 

1,0*5 

throttle 

•ETTIN* 

TO 79,0 

PERCENT IN 5 

<« 


71,0 

•coin 

11514,0 0,0 770,0 

TURN TO HEADING 45,0 DEG, 

162,5 

1*0,0 

,204 

,01 

1.010 

7*,0 

11744,4 

.0 

*14,* 

162,3 

140,1 

.204 

.01 

1,007 

71,0 

1404(,1 

1,2 

444,0 

162,3 

140,1 

,204 

.01 

1.001 

74,0 

14140.9 

3,6 

**0.1 

162,3 

1*0,* 

,204 

,00 

1.015 

79,0 

14412,4 

16,7 

*00, S 

162,3 

140,1 

,204 

.01 

1,010 

74,0 

14009.2 

56,2 

412,8 

162,3 

1*0.1 

,204 

.10 

1,041 

77,0 

19194,4 

67,6 

494,4 

162,6 

1*0.1 

,204 

.11 

1,172 

74,0 

194*4,2 

114,3 

474,4 

162,6 

1*0.1 

.204 

.12 

1,221 

74,0 

19442,1 

1TG # 2 

4*2,7 

162,5 

140,1 

.*«7 

.** 

1,2*1 

•0,0 

19494,2 

259,2 

10*4,1 

161,6 

140,0 

.207 


1,*14 

•1,0 

1*110,4 

557,1 

10l*,4 

163,1 

1*0,7 

,207 

»*9 

1,219 

•*,0 

1*441*4 

471,6 

1018.7 

163,5 

1*1,1 

,20* 

,« 

1,207 

•i!o 

14704;* 

602,5 

10*1,1 

164,1 

141,4 

,204 

1,01 

1,144 

•4,0 

14*40,4 

766,6 

10**, 0 

164,6 

1*2.2 

,290 

1.2* 

1.141 

•9,0 

171*7,4 

910,5 

1071,4 

165,3 

1*1,4 

,*91 

1.31 

1,1*1 

•4,0 

i7i«*;i 

1044,4 

10*1,0 

166,3 

1*1,7 

.*92 

».w 

1,172 

07.0 

1794*; * 

1269,3 

10*1,4 

167,6 

uo.o 

,291 

1,10 

1,1*4 

••,0 

17109,4 

1461,2 

1104,2 

167,6 

140,4 

,29* 

,*0 

1,15* 

S6 t 0 

1*004,2 

149*, 1 

m *;« 

166,0 

1*9.1 

.299 

.S» 

1,14* 

90,0 

i* 2 io;* 

1*97,4 

!1ST,T 

166,2 

1*9,0 

,299 

.1* 

I.14B 

91,0 

1*411,0 

2097,4 

1194,2 

166,2 

149,0 

,299 

.01 

1,101 

91,0 

l44lli* 

2297,4 

1 l*t,4 

166,2 

1*9,1 

,* s * 

,00 

1,044 

91,0 

1*411,1 

249*, 2 

1204,* 

166,2 

149,1 

,*95 

.01 

1,050 

94,0 

14011 9 

249*, 4 

1211,1 

166,2 

149,* 

.299 

.01 

1,044 

95,0 

1*211,4 

14*11,7 

2*94,4 

1199,6 

166,2 

149,* 

,295 

.01 

1,044 

96,0 

109*, 4 

12*0,1 

166,2 

149.1 

.254 

.01 

1,046 


CD ALPHA 

• ANNA 

R/C 

LOAD 

THRU6T 

FU6 § 

ROLL 

HZ AD HQ 

(DEG) 

(DEO) 

CPPN) 

FACT 

(LSI) 

ANG' 

ANGLE 

(DEG) 


'1363 

a, 10 

1,10 

401,1 

1,10 

37126 

10,1 

0,0 

0,0 

,1267 

8,00 

1,«* 

1101,2 

1,10 

37106 

10,4 

0,0 

0,0 

,1232 

7,40 

4,9* 

1300,0 

1.10 

37)66 

11,5 

0,0 

0,0 

,1166 

7,00 

9, *7 

1900,4 

1.10 

37)76 

12,1 

0,0 

0,0 

.1105 

7.79 

9,*7 

1702,1 

1.10 

37)64 

U.T 

0,0 

0.0 

V 1 Q 62 

7,79 

4,47 

1*02,4 

1.10 

37)36 

n .« 

0,0 

0,0 

'1016 

4,40 

7,1* 

2044,4 

1.01 

37)17 

19.1 

0,0 

0,0 

,0667 

4,90 

7,15 

2049,5 

1,00 

37)16 

!*.• 

0,0 

0.0 

,0676 

4.50 

7,1* 

210*,* 

1,00 

37)36 

t *,4 

0,0 

0.0 

'0621 

7,50 

7,41 

2114,4 

1,00 

37)16 

19,* 

0,0 

0,0 

v 0670 

4,90 

7,47 

2111,5 

.’4 

37)30 

15,0 

0,0 

0,0 

'0611 

4,90 

7,44 

2124, * 

,4* 

37)66 

15.* 

o.o 

0,0 

'0675 

10,90 

7.91 

2144,4 

1,04 

37)36 

17,0 

0,0 

0,0 

'0662 

10,19 

7, ** 

2*47,4 

1,01 

37)37 

17,* 

0.0 

0,0 

,0637 

10,05 

8,00 

2247,4 

1,00 

37136 

17,0 

0,0 

0,0 

'0616 

4.49 

8,04 

2100,7 

.♦4 

57131 

17,0 

0,0 

0,0 

.0626 

4,49 

*,04 

2304,0 

,49 

17134 

17,0 

0.0 

0,0 

,0616 

4,05 

4,07 

2104,4 

.0* 

37)14 

17,0 

0.0 

0,0 

,0626 

4,45 

8,07 

2104,4 

.04 

37)33 

17,0 

0.0 

0,0 

t 0 B «6 

4,45 

*,07 

2307,1 

,0* 

37)32 

17,0 

0,0 

0,0 

t 0626 

4,45 

8,04 

2305,7 

.** 

37)3! 

17,0 

0.0 

0,0 

,0626 

4,45 

8,04 

2105,4 

.4* 

37)30 

17,0 

0,0 

o;o 

,0775 

4,10 

7,41 

2247,4 

,04 

33630 

1*.* 

0,0 

0,0 

,0766 

0,45 

7.47 

2114,7 

.01 

33673 

19.4 

a .0 

,0 

.0766 

0,40 

4 , ** 

19*0,0 

.♦0 

32)16 

14,* 

7.0 

.9 

,0755 

4,05 

4,11 

17*0,1 

.01 

30136 

>«.» 

12,0 

l.» 

,0771 

4,29 

5,50 

1574,5 

.02 

16402 

19, • 

17.0 

9.1 

,0637 

10,09 

4,81 

117*,* 

,0* 

27644 

19,* 

*2,0 

9,2 

>15 

11,05 

4,32 

1200.3 

1.09 

27643 

14.4 

27,0 

•ft 

,04*5 

11,70 

1,4* 

1142,7 

1.10 

2764) 

1«.7 

90,0 

11,4 

,04*9 

11,70 

1,47 

1055,1 

1,10 

27636 

1*.4 

90,0 

19,9 

,0401 

11,45 

3,35 

945,1 

1.10 

27633 

1 M 

90,0 

1*,0 

,0474 

11,40 

1,04 

475,0 

1.10 

27626 

19,4 

90,0 

22,7 

,0440 

11.90 

2,72 

7*4,1 

1.10 

276)6 

19,1 

90,0 

**.« 

,0497 

11,40 

2,40 

444,9 

1.10 

*7604 

II .* 

90,0 

90,1 

,044* 

11.10 

2,0* 

*04,1 

1.10 

27766 

I*,* 

90,0 

99,7 

,0410 

11.20 

1,41 

999,7 

1,10 

*777* 

12,1 

19,0 

»M 

,0425 

11,09 

*.01 

540,4 

1,10 

*7714 

t*.l 

20,0 

94,0 

,0*14 

10.49 

2,3* 

4*9,4 

1.10 

*7736 

l *, 9 

1».0 

• t ,0 

,0404 

10,89 

2,80 

*2* ,7 

1,10 

177*4 

>*.‘ 

10,0 

09,4 

,0*4* 

10,79 

1,3* 

1009,1 

1,10 

177)4 

19,1 

»;• 

44,9 

,0*44 

10,79 

4,02 

1145,4 

1.10 

*7706 

19,0 

*.0 

•4,0 

,0445 

10.19 

4.9* 

1144,4 

1,05 

1770* 

19,7 

0,0 

#»,0 

,0407 

4,70 

4, *3 

1414,1 

1.02 

17707 

19,9 

0,0 

49,0 

,074 j 

4,90 

4,41 

145#, 0 

1,00 

27707 

19.4 

0.0 

49, P 

,07*7 

4,09 

4,41 

1444,9 

1,00 

27T0T 

19.4 

0,0 

«».• 

,07*7 

4,*9 

4,44 

1444,4 

1,00 

27706 

I *.* 

0,0 

41,0 

,07*7 

*.«» 

4.49 

1470,1 

1,00 

27706 

19,4 

0.0 

•9.0 



TAKEOFF CONTINUED 


TIHE 

X DI8T 

V D 1ST 

ALT, 

TA8 

EAS 

MACH 

ACCEL 

CL 

t«c> 

(FEET) 

(FEET) 

(FEET) 

( KTS) 

<*T8) 

NO. 

(f PS2) 


97,0 

19611, « 

1 

1304,9 

168,3 

163,1 

,256 

.01 

1,049 

98,0 

19*12,0 

5459,4 

1329,1 

168,3 

163,0 

,256 

.01 

1.046 

99,0 

20012,1 

5659,6 

1353,9 

168,3 

163.0 

,256 

,00 

1,050 

(00,0 

20212,2 

1659,9 

1378,1 

168,3 

16<I,9 

,256 

,00 

1,050 

101,0 

20412,0 

4060,2 

1402,6 

168,1 

164,8 

,256 

,00 

1,050 

101,0 

20912,9 

4260,4 

1427, 1 

168 c 3 

164,8 

,256 

,00 

1,050 

(0»,0 

20*12,7 

4460,7 

1451,6 

168,3 

164,? 

,256 

,00 

1,050 

(0«,0 

21012,9 

4661 , 0 

1476.1 

168,3 

164,7 

,256 

.01 

1,050 

109,0 

21215,0 

4861,3 

1500,5 

168,3 

164,6 

,256 

• 01 

1,050 

109,0 

21015,2 

5061,5 

1525.0 

168,3 

164,6 

,256 

.01 

1,050 

107,0 

21913,4 

5261,8 

1549,3 

168,3 

164,5 

,256 

,00 

1,050 

100,0 

21»13, 9 

5462,1 

1573.9 

168,3 

164,4 

,256 

.00 

1,054 

109,0 

22013,* 

5662,4 

159* ,4 

168,3 

164,4 

,256 

.01 

1,050 

110,0 

22214,0 

5862,7 

1622,9 

168,3 

164,3 

,256 

.01 

1,050 

111,0 

22414,2 

6063,1 

1647,3 

168,3 

164,3 

.256 

,00 

1.058 

112,0 

22614,4 

6263.4 

1*71,7 

168,3 

164,2 

• ?56 

,00 

1,05* 

1 1 5,0 

22*14,6 

6463,7 

1696,2 

168,3 

164,2 

,256 

,00 

1,05* 

n«.o 

23010,0 

6664,0 

1720,6 

168,3 

164,1 

,256 

,00 

1,058 

115.0 

23215,0 

6864,4 

1745,1 

168.3 

164.1 

, 256 

.01 

1,05* 

ADVANCE THROTTLE 

setting 

TO 95, 

0 percent in 

3.3 SEC, 


11*. 0 

23415,2 

7064,7 

1769,* 

198,3 

164,0 

.25* 

,04 

1,113 

1 1 7, 0 

23615,4 

7264,9 

1795,* 

196,3 

164,0 

,256 

,01 

1.156 

lt«.0 

23*15,5 

7465,0 

1 * 25 ; 1 

198,4 

163,9 

,256 

.04 

1,190 

ii’.o 

20015,0 

7664,8 

1*57,6 

198,4 

163,8 

,256 

.01 

1,139 

120,0 

20210,5 

7864,4 

1*92,* 

19*. 4 

163,8 

,256 

,01 

1,0*1 

121,0 

20015,9 

8063,9 

1929,0 

198,4 

163,7 

,256 

.01 

1,092 

122,0 

20*13,3 

8263,4 

19*5,7 

198,4 

163,6 

,256 

.01 

1,050 

125.0 

24*12,6 

8462,9 

2002', 4 
OF 250 
2039,0 

198,4 

163,3 

,256 

♦ 01 

1,050 

ACCELERATE TO CLIKB SPEED 
120.0 25012,0 06*2,0 

.0 

198,5 

163,5 

,256 

,20 

t.ois 

125,0 

25211,7 

*** 2,2 

2074,5 

198,7 

163.6 

,257 

.51 

,975 

12‘,0 

25411,9 

90*2,5 

2107,7 

199,1 

163,9 

,257 

,95 

,919 

127,0 

25915,1 

9263,* 

2137,4 

199,* 

164,5 

,259 

1,49 

.897 

12», 0 

25*15,6 

9466,4 

2192,9 

170,9 

165,5 

,240 

2,07 

,873 

129,0 

26019,9 

9670,* 

21*3,5 

1T2;3 

166,8 

. 292 

2,66 

,850 

150,0 

26229,2 

9*77,3 

2199, 0 

174,0 

168,3 

,295 

3,26 

,830 

151,0 

26455,1 

100*9,5 

2209,3 

179,0 

170,4 

,29* 

3.24 

,976 

152,0 

26646,3 

10297,9 

221*, 1 

178,0 

172,2 

,271 

3,25 

,956 

135,0 

26*59,9 

10511,3 

1227,0 

179,9 

174,1 

,274 

3.29 

,917 

150,0 

27075,7 

10727,3 

2235,8 

1*1,* 

175,9 

,277 

3,28 


155,0 

27295,9 

10945,6 

2244,6 

1*3,* 

177.8 

,2*0 

3.29 

,900 

begin 

139,0 

TURN TO HE*DING -15 
27514,4 111*9,2 

.0 DEG, 
2253,5 

1*5,7 

179,7 

,2*3 

3.29 

,8*1 

157.0 

27737,9 

11388,0 

2292,3 

1*7,7 

181,5 

,289 

3,29 

,871 

15G.0 

279*6,3 

11910,3 

2271.1 

1*9,6 

183,4 

,2*9 

3,29 

,866 

159,0 

2*201,6 

11*29,7 

22*0,0 

191,5 

1«5,2 

,292 

3.20 

.867 

100,0 

2*445,7 

12044,2 

228*. * 

193.4 

187.0 

,29(1 

3.12 

.877 

141,0 

2*700,3 

12251,1 

22*7,6 

195,2 

188,7 

,297 

3,09 

,8*0 

102,0 

2*966,7 

12447,6 

2306,1 

197,1 

190.5 

,300 

3,12 

,865 

103.0 

29244,4 

12933,* 

2314,4 

198.9 

192,2 

,303 

3,12 

,852 

100,0 

29532,5 

12*09,1 

2322,7 

200,7 

194,0 

,309 

3,12 

.837 





END OF TAKEOFF 


Cl 

CD 

ALPHA 

GAMMA 

R/C 


(DEG] 

(DEG) 

(fPH) 

,B22 

,0565 

6,60 

1.5* 

499 1 6 

,804 

,0571 

6,57 

l.J» 

500 . 1 

f 795 

,0564 

6.25 

1.37 

500.1 

,778 

,0532 

6.75 

1.5* 

500*0 

,764 

,0503 

7,2* 

1.55 

500*1 

,751 

,0492 

7.U 

1,33 

900*3 

,715 

,0460 

6,6 1 

1,32 

500.4 

,723 

,0472 

6,7* 

1.51 

500*1 

,709 

t 0437 

6,72 

1.50 

499,3 

,896 

,0425 

7,37 

1.2* 

500*1 

,685 

,0416 

7,22 

1.27 

500*3 

,672 

,0407 

7,06 

1,2* 

500*5 

,660 

.0398 

6,90 

1.25 

499,9 

,669 

,0390 

6,77 

1,24 

500.4 

,6J9 

,0363 

6,63 

1,23 

499,6 

,626 

.M5 

,0375 

;0366 

t;tt 

1:1! 

ItS;} 

, 59 1 

,0151 

6,03 

1,20 

500,4 

,565 

,0314 

5.69 

1.1* 

500,4 

,546 

,0121 

5,92 

i.i; 

500,5 


,0112 

5,22 

1.17 

500,4 

.916 

,0105 

5,07 

1,16 

500,7 

,507 

,0101 

4,96 

1.15 

500,7 

, 499 

,0296 

4. *6 

1,14 

500,5 

,509 

,0102 

4.9* 

1,20 

534,9 

,519 

,0107 

s.n 

1,43 

640,7 

,529 

,0311 

5.23 

l,*l 

*1*,0 

,519 

,0118 

5.16 

2,34 

1064,8 

,546 

,0124 

5,48 

3,02 

1376.5 

,554 

,0326 

5,56 

3, *2 

1752,6 

,552 

,0127 

5,53 

4,66 

2143,4 

,550 

,0325 

5.51 

5,50 

2533,4 

,346 

,0124 

5,4* 

6,34 

2921,5 

,511 

,0101 

5.01 

7,11 

3276,5 

,470 

,0261 

4,4* 

7,38 

3401,1 

,460 

,0276 

4,36 

7,45 

3432,9 

,436 

,0273 

4,31 

7,47 

3442,1 

,497 

,0274 

4,12 

7,47 

3443,6 


LOAD 

thrust 

FUS, 

ROLL 

HEADNG 

PACT 

(1881 

ANG, 

ANGLE 

( DEG) 

1.10 

14097 

7.0 

•24*3 

27,6 

1*09 

34044 

8.8 

•24*1 

25.4 

i.io 

11990 

6,6 

•24*7 

22,7 

1,10 

11917 

7.1 

•24,2 

20,5 

1.10 

11861 

7.6 

•24*7 

1*,0 

1.10 

11629 

7,® 

•24,6 

15.6 

1.09 

13769 

7.4 

•25,0 

15. » 

1.10 

31712 

7.1 

•24*7 

10,0 

1,09 

33656 

7,0 

• 21*9 

*.6 

1.10 

13596 

7.7 

•25*2 

8.* 

1.10 

13539 

7.5 

•24*6 

o.l 

1.10 

13460 

7.5 

•24*5 


1*09 

13422 

7.1 

-24*0 

■«5 

1,10 

31161 

7.0 

• 24 f 5 

-4,4 

liio 

13104 

*.* 

• 24*3 

■4 * 6 


llffl 

M 

■ 24*6 
• 24t0 

:W 

1,07 

11126 

*.4 

•21.0 

• 10,0 

1*04 

11066 

5.* 

• 16,0 

•12,4 

1.02 

13004 

5.* 

-11*0 

-15,5 

1,00 

32942 

5,0 

•6*0 

• 14,2 

1.00 

12660 

5.4 

• 2*0 

• 14,5 

1.00 

12816 

5.1 

• 2.0 

• 14,6 

1,00 

32756 

5.0 

•2.0 

• 14,8 

!«01 

32695 

».4 

•2,0 

-1S.0 

1,07 

32616 

5.5 

0,0 

-15,0 

1*11 

32562 

6.0 

0.0 

• 15,0 

1.19 

32511 

‘.7 

0.0 

• 15,0 

1,17 

32490 

7,5 

0,0 

-15,0 

1*20 

12456 

8.4 

0,0 

-15,0 

1,20 

32429 

*.4 

0,0 

• 15,0 

1,20 

32411 

10,0 

0.0 

• 15.0 

1.20 

32400 

10,* 

0,0 

• 15,0 

1.11 

32397 

11,1 

0,0 

• 18,0 

1.02 

32595 

10,0 

0,0 

• 15,0 

1,00 

32194 

10,* 

0,0 

»is,o 

,99 

12192 

to,* 

0,0 

• 18,0 

.99 

32191 

10,* 

0,0 

• 15,0 



APPENDIX C - PROGRAM LISTING 


T AK0 F F 

SUBROUTINE TaKOFF flNPC , TDCN, HGROSS , S*I NG # XENG# VR# VEND ) 

C 

C SUBROUTINE TAKoFF COMPUTES THF TAKEOFF maneuver of a given aircraft# 
c INCLUDING GROUND ROLL AND CLIMBOUT, 

£ provisions are made for changes in flap# vectored thrust angle and poker 

C SETTINGS AS FUNCTIONS OF SPEED AND ALTITUDE* 

C 

C HEADING ANGi.ES ArE DETERMINED RY ALTITUDE AND GROUND DISTANCES, 

C 

C FOILING $0 H E COMMENT CApOs Will. BE TWO NUMBERS IN p ARE Ny HE S I $ WHICH 

C will give the approximate statfment^numbe* RANGE OF THE FUNCTION 


DESCRIBED IN ThAT COMMENT, 





EXTERNAL DERIVi 
EXTERNAL DKRIV2 
REAL M(J ( M tlTfRXC 
COMMON /UMJV/ npc 

t NSC 

,IDC 

t h #sr f r 

#w 

r 

twF #Fm #VMO 

, EmmO 

» ALPHLO 

#CLALPH#SN #ar 

t B 

t 

<?E YE W ,£nP f T A 

3# KS I ZE 

t WG 

, wGS 

#KnRITE#DLMC4 




COMMON /AFRO/ VEL#0S/HA8S#THRUST # TVECT# angle# delfd#delspl# alpha# 
9CVfCY#CL#CD#RH0#DRC0#IFAST 
COMMON /POLL/ PHIfttOCMIN 
COMMON /EXC-HNG/ SRC1LL# 93*5# y3B#T5j 
COMMON /XPnLL/T<30) ,Nf Q#MU# NREV 

common /XFLATE/ sf /5),R0C#THFMAX#THETAF#XLF#)(LFMAX#K 
COMMON /UMT/IUNJT.G 
DIMENSION ans(A1 

DIMENSION XDELKOC5) #XHFLAP(snxVFL AP(S: 

DIMENSION XHRHRC5) #yvPhR( 5) #XP0WfcR(5) 

DIMENSION YNUfj' 1 #XHVE'CT<5! #XVVECT(5) 

DIMENSION XPANGE(S), XHhEAD(B) # XHEADfS) 

YYYtXXX# Yi#>2* Xl#X?T » Yl f ( Y2- Y 1 ) * ( XX X-X 1) / ( X?-X X) 

DATA xDELFU/j5,0 f S. 0,2,0, 0*0, 0,0/, XHfL a P/O , 0 , *5 0 , # 0 , 0 # 0 , 0 , 0 f 0 / # 
^XVFLAP/O* 0*0,0,200, ,210*, 0,0/, XPO^ER/1 ,0« 1,0# l ,0# 1 ,0#l,0/, XHpWR/ 
9Q,0, 0*0,0, 0,0,0#0»0/# XVPwa/Q,0#999,,999 fJ 999,#999 t /# XNU/0,0#0 # 0# 
9D,0#0*0#0 B 0/, XHVECT/O.0 f 0,0,0 ( 0#0,0#0,0/# XVVECT/0,0, 999 , # 999 , # 
9999, #999, / 

DATA XPANGE/5* 100,00/# XNhf AD/5*99999, /, XHE A D/5*0,/ 

DATA CDGEAR#nKLPl>T#DTABS#DTGR#OtPDWN#DTPJP#DTVECT#EYEMNG#MAPT# 
9HGR#HMAN# TOUT #UM#PmaRG#RTCL# THTFLV# THTSCP/O,0#3,0#0 t 0#s # 0#5,0#6 t 0# 
<M o, #1*0*0, 0,25, 0# 1000, #0, 0. 02, 0,00# 75 0,, 15 ,,10,0/ 

D a TA ROLLMX#RnLRAT,RnCMTN#HDT/l 5 *0,5,0 F 250, #35,/ 

DATA METRTC/IHm/ 

NAMELTST/NAMi/CDGF.AR# DAOy#DF LPOT# OT ABS# DTGR# OTPDnN# DT PUP# DI VECT # 
9EYFWNG#HAPT#hGW#mman,HhAX# I OUT , UM, NP AGE# PMARG • RTCl , TNTFLY# 

9THTSCP* XlFmAX#HOllmX, RUI.RAT, ROCMIN/HDT 
NAMFLIST/NAM^/x OF LED#* HF LAP, X V E L A P , X P D itf E R # XHPwR, X VP WR # XNu# XH VE C T # 
9XVVECT 

NAMEL IST/NAM3/XRANGF #XHHEAD# XHEAD 
GO TOCS, 6, 5) , INPC 


36 

5 nPAGE 

* 33 

37 

dadt = 

M 

40 

hmax s 

50no, 

42 

XLf MAX 

= J.IO 



T AKOFF 


C 

c 

c 


SFt -jp f 0® IJ*T7J> Cn-gVF SSlON 


04 

04 



, AND / N AM3/ 



TAKOFF 

166 IRQU » 1 

166 MROLL ■ 0 

167 IHEAO a t 

170 PHIMAX » R0L.LMX/5T ,3 

172 TINT ■ 10,0 

C 

C find gear DRAG If DEFAULT OF CDGEARaO.O USED 
C 

174 IFCCOCEAr ,Eq. 0,0)CDG£Ar ■ ( 0 , 0032/SW)* W**0 ,80 

204 CDGEAR a CDGFTR*COGEAR 

206 URCD a COGEAr 

206 IF ( IDCN ,NE. 9 ) GO TP 249 

210 TEMP » 54,0 ♦ PTABS 

212 MG a WGROSS 

213 WOS • W/$w 
C 

C FIND STATIC THRUST /WE T GHT 
C 

215 CALL EnG I nF 10,0,0.0,0,0, 1,0, wF»KEnG) 

221 TON a ENP * THRUST / W 

C 

c write out program inputs 

c 

224 WRITE(6,20I)HAPT,TEMP . _ 

201 FORMAT!//, 36H ** INPuiS TO TAKE OFF - ALTITUDE « , F6 , 1 , 2X , 1 4HTEMP 
4ERATURE » ,F5.1,7H DEG, ) 

243 WRITE(6,2021WG,SW, THRUST 

202 FORMAT!/, 6X, J4HA/C CHARACTERISTICS, /,4X,17HGROSS RAMP WT, a ,F8,0, 
43 x,12HwING Area a , Fe , 0 , 3 X , 26MST ATlC SE* LEVEL THRUST * ,F6,0> 

260 WfllTE(6,203)WDS,TOW 

203 F0RMAT(Rx,l5HWlNG LOADING a ,F5,I ,3X, 16HTHRUST/WEIGHT a ,F4,3) 

273 WRITE(6,204TENP, CDGEAR, EYEW,THTSCP 

204 FORMAT!/, 6X, 15HA/C PARA ME TF RS, / , 4X, 1 4HN0, ENGINES * ,F3,l»J*r 4MC 
9DGEAR a ,F6.4, iX,9HETEWNG a , F4 , 1 , 3X , 20 MT A II SCRAPE ANGLE « ,F4,|) 

312 WRlTE!6,205>XLFM*x,HGR,THTFlY,HMAN,RTCL 

205 FORMAT !/,6X,30HFL JGHT PATH CONTROL PARAME TERS, / , 4X , 1 8HMA X LOAD FAC 
9TOR a ,F4.2,3X,?3HGEAR RETRACTION ALT, a , F5 , 1 , 3X, 18HMAX FLOOR ANG 
4LE a ,F«. 1,/,9X, 16HMANEUVER ALT, a , F5, 0 , 5X , 27HACCELERATE RATE OF 
9CLIMB a ,F5.0) 

333 WRITE! 6,206) DADT, DFLPDT,DTGR,0TP0WN,DTPUP,DTVECT 

206 FORMAT!/, 6X.25HPARAMETER VARIATION RATES, /,4X, 7HDADT a ,F4,l,3X 

9 , 4HDF LPDT j , Fa, | , 3X, 7H0TGR s , F 4 , 1 , JX , 9HOTPDWN a , F« , 1 , / , 9X, SMOTO 
9UP » ,F4,1 ,3X,9H0TVECT a ,F4,1) 

356 WRITE(6,207) 

207 FORMAT!//, fcX,42MP0WFR, VECTORED THRUST, ANO FLAP SCHEDULES) 

365 WR1TE(6,20B)1XP0WER(I),I ■ 1,5) 

208 FORMAT!/, 9V,22HTHR0TTL£/P0WFR SETTI NG , / , 12 X , 8HPWRSE T ,5F9,2) 

376 NRITE(6,209) (XVPWR(I),I* 1,5) 

209 F0RMAT!1?X»8HSPEE0 ,5F9,1) 

407 WRITE!6,210)(XMPWR(I),I * 1,5) 

210 FORMAT(12X,RHALTITUDE,5F9,0) 

420 WRITE(6,2tl)(XNU!l),l a 1,5) 

211 FORMAT!/, 9X,2|HVECT0RFD THRUST A NGLE , /, 1 2X , 8H ANGLE ,5F9,1) 

431 WRITE(6,209)(XVVECT(I),I ° 1,5) 

436 WRlTE(6,2lO) CXHveCT!!), I » 1,5) 

447 WRlTE!6,2l2)!XDELFD!I),I * 1,5) 



TAKOFf 


SF9 # 1) 


L n 


46 a 

*71 

50? 


5)5 

524 

535 

546 

557 

557 


572 

6)0 


626 

627 

631 

632 
634 

634 

635 

636 

636 

637 
640 
642 


642 

644 

645 


2 i 2 FORMAT C/,9X,2tHfLAP DFFl.EC T I ON A NGLE # / , 1 2*, 8HDELFD , 

WRITEC6#209) CXVFtAPf Ii# I * 1#5) 

WRITE(6,210) (XHFLAPiI), l « 1,5) 

213 F0RMAT(A9x!76HALt SPEEDS ARE INDICATED AIR SPEEDS AND ALL ALT I TUD 
9ES ARE AB801UTE ALTITUDES) 

wRI TE 16# 2 1 4 ) 

214 F0RNATC//,6X r 18NDFPARTURE HEADINGS#/) 

W«ITEC6#2)5) (XRANGF (I) ,7*1 ,5) 

2)5 FORMAT C 12X,8mRANGF ,5FlO f l) 
wR!TE(6#2)6) (XHHEAOd),! 3 1,5) 

216 FORMAT (12X#8H ALTITUDE ,5F10,D 
WRITE (6# 21 7) (XHEAncn , 1*1 ,5) 

2)7 FORMAT (12X#RHH FADING ,5F!0 a |) 

299 CONTINUE 

IFdOCN ,FQ, 9) 

4WRITE (6,999)HAPT 

999 FORMAT ( )H 1 9 //Z 1H TAKEOFF (flfVATION *,F6 t 0#4H FT)#//) 

IFCIOCN .EO. 9 .AND. IUNIT f E0. 1) MR I Tt C 6 # 1 0 0 0 ) 

IF { IDCN .FQ. 9 , A NO t IUNIT t en„ 2)NRlTEC6 i 2000) 

l SET FLAPS# vECfORED THrUST angle and POWER SETTING FOR GROUND run 

pWHSET * XpnwERfl) 

ANGLE * XMUU) 

DFLFO * xntLFDO) 
t hfmax * thtscp 
MABS * 0.0 
vFL * 0, 

ZFRO * 0, 

7 * 0,0 
THFTAF * 0. 

QS * 0,1 
S(7) * HAPT 

e« * Of 

C GROUND POLL INTEGRATION vARIA0L.ES 

C TCI) * NUMBER OF equations 

c T ( 2 ) ■ TIME (SEC.) 

c T ( 3 ) 3 T IMF INTERVAL, STEP SIZE (SEC.) 

C T ( 4 ) = VFICCITV (FT. /SEC.) OR (M/SEC,) 

c t< 5 ) * distance ift.) or (m) 

c T(6) * ACCtlEPATTON (FT f /SEC t ** 2 ) OR (h/SEC.**2) 

C 

TUT * 2 

T ( 2 ) = 0.0 

T 1 3) « «•) 


646 

646 

c 

IS?! : 8: 


C 

C 

obtain atmospheric variables 

650 

c 

CA|_L ATMO$(HAPT#OTARs,AnS) 

652 


SA s ANS<4) 

653 


PHO = ANSI 3) 

655 


NCOIINT 3 0 



TAKOFK 

655 


ALPHA = EYE* 


04 

O' 


657 

665 

666 
667 
701 
703 


767 

77 1 

772 


773 


777 

1011 


1012 

1010 

1015 

1017 

1023 


1047 

1052 

loss 

1056 

1056 

1060 

1062 

1066 

1073 

1074 

1075 
lino 


c 

c find Initial conditions uf thrust and force coeff, for ground roll 
c 

CA^l EwGjnF (S( 7 ),DTA&S,Fm,PKRSET, WF#*EnG) 

CALL AROOYM 

CALL INTS(T,2#2, 1., U, Ur 1., t ,, 1,# DER1V1) 

ta s thrust * enp 

IFIIOCN «FQ* 4) 

♦WRITE (6 , 1002)1(2 )! T f 5 ) , ZERO, ZERO , ZERO , ZERO , ZERO, T ( 6 ) , CL, C D f AL PH A , 
92ER0 # ZERO, ZERO# TA,*FRO, ZCROi ZERO 
C 

c main ground roll Integration Loop <im 20) 
c 

l n Count c ncount f 1 
ai.phaj v alpha 
thetaj s thftaf 
c 

C CHECK IF LIFT s height - LIFTOFF 
C 

IF(03*CY .GF. «OGO TO 120 
C 
C 

C CHECK SPEFD FOR YR 
C 

106 IF ( T ( 4 ) * vELFTR*SORT(AnS (7 J ) . GF. # VR)GO TO 101 
GO TO HO 
C 

C begin rotation 

c 

101 IFCJJ2 .EO, 2)Gn TO 103 
JJ2 = 2 

yK T S 3 T { T * VELFTR 
E AS = YK1$*S0PTCANS(71 ) 

IF(TDCN # FG. R)«miE(b# 1009)T(2), VKT$,EAS 
t 00 R FORMAT ( 1 X, 1 RMRUTATIflN {TIMf a , f 6 , 1 , 3X , 1 0 H A NO TAS b #F 7,1,2X, 
96HEAS * # F 7 • 1 , 1 H ) ) 

C 

c increase alpha at rate oadt and restrict to value ,le, tail scrape 
C angle ( 103-102) 
c 

103 Alpha s alpha ♦ tc 3 )*dadt 

IF ( (ALPHA - FVEw) , LE , lHEMAx)GO yO 102 
ALPHA a THEHAX ♦ FYFW 
GO TO 1!0 

102 THfcTAF a ALPHA - FVF* 

110 F« s T ( 4 ) / S A 

OS = 0.5 *RhO*Sw*T ( 4) *T(<J) 

CALI ENGINE (S( n ,DTA0 S,EM,PWRSET ,wf ,KENG) 

1 F AST = 0 

VEL = T ( 4 ) 

TA = EnPaTmRuST 

CALL iNTMfT,?, 2, #!.#!•#! 


C 

C 


HAKE integration stfp 


OERIVIT 





TAKOFF 


1111 

1115 

1124 

1133 

1H5 


1153 

11)6 

llan 

1143 

1143 

114b 

1152 

1160 


1237 

124t 

1243 

1243 


1254 

1272 

1310 

1311 
1313 
1317 
1321 
1350 


W * w - nF^EN^TfD/SbOO, 

If(T(2) «LE. 40.1G0 TO 113 
kRITt(6,1057) 

1057 fORMAT(/,35H TIME LIMIT FOR GROUND RUN EXCEEDED) 
RETURN 
113 CONTINUE 

: COMPUTE ALPHA oot and theta ddt 


y DAD T S {ALPHA - AL PHA J ) / T ( 3 ) 

UTHTOT = (THETAF - THET A j) /T ( 3) 

IF(NCOuNT ♦LT. 10)GO TO 1 
NCOUNT * o 

VKTS * T < AM * VFtFTR 
E AS = VKTS * SQRT(ANS(7)) 

JFtTDCN ,NF,9)&0 TO i 

♦ iRITE(6 r lOO?)TC2)fTC5) # 2#ZERO,VKTS,rEAS f EH # T ( 6 ) , CL , CD , ALPH A , ZERO t 
^ZFROf 7fR0 # T A# THETAF, ZERO, ZERO 
I PAGE = IPAGfc +1 
TFUPAGE .LT . NPAGF1GO TO 1 
IPAGE = 0 
JFUOCN , FQ # 9) 

♦WRITE (b# 99b) 

99ft F0RHAT(1W1 , /, 19H TAKEOFF C QNT 1 NU t D • /) 

JF { TDCN .FQ, P ,AND. [UNIT , EQ , 1 ) MRI T£ ( 6, 1 000 > 

1 F ( I DC N ,EG # 9 .AND. I UNIT # EQ, 2)*RITE(b,20Q0) 

GO TO 1 

120 V'KTO s T(O) * velftr 

FASTO * VKTH * SQRT(ANS(7)) 

SROLL * T (5) 

JF ( IOUT # E0g 1 ) EnP = EnPOUT 

IF ( I DCN ,FU, 91 WRT TF. ( b, 1 0 1 0 3 T ( 2) , T (5) , VKTO f EASTO 
1010 FORHA T ( l y , 1 THE I FTOFF A J 1 H€ = » Eb . 1 , 2X # 7HD I ST s f FB f U2X 
9,6HTAS = ,FM,1X,*HMS = , *7,1,1*)) 

1000 F0RMAT(13|H UME X DIST Y DIST ALT f ntrl T *f H Pulj S FUS RO 
9CCFL CL CD ALPHA GA*MA R/C LOAD THRUST *US, RO 

91 L HEADNG,/, 131H (SEC) (FEET) <KfET) ( p tET) ( KTS ) CKTS) 
9NO f (FPS2) (DEG) (DEG) (FPM) FACT (L93) 

9. ANGLE (DEG),/) 

2000 FORMAT ( 1 31 H TIME X DIST V DIST *LT« f Rn 

■ — * gamma r/c load thrust rus, hd 


ang 


MACH 


9CCEL 


CL 


to 


alpha 


(MTRS) (MTRS) (MTRS) (M/S) (M/S) 
(DEG) (DEG) (M/M) FACT (NTS) ANG 


9Ll HE AONG, /# 131 h (SEC) 

9nQ t (MPS2) 

l0fi2 , rORH»T(t x ,FS^WF9^1,F<»,l,F‘>.1^7,l,F6,l,F6,J,P7 1 2,F7 t 3,F7,«f2F7.2, 

9Fe,l,F6,2,F9,0,F6.1.F 7.I,FU,1) 

: fcNO of CROWD ROU. - PtRJN airborm* portion OF TAKtpFF 


1JSJ 

ISS5 


VMARG s PMAPG*VFnD 
THfcMAX a THTFLV 


FLAP# Angle and power schedules set up FqR airborne PORTION 



TAKOFF 

1356 HflAp s XHHAP(2) 

1560 VFLAP s XVFLAP(2) 

1361 HPOWF* * vhP*R<2) 

1363 VPOWF.R = XVP*R<2) 

1364 HVECT = XHVFC.T (21 

1366 VVFCT a XVVl:'CT(2) 

1367 RANGE s X R A NG F { 1 ) * R M G F T R 

1371 HHEA[) * XHHFAO(l). 

1373 PHI = 0,0 

1373 XlF a 1.0 

C 

C AThHoRnL InTERGPATION variarlf$ 

C S(l) * NUMBER OF fouations 

C S(2) * TIME (SEC,) 

C S ( 3 ) * time; INTERVAL, STEP SIZE (SEC-) 

C S ( 4 ) = VELOCITY along FLIGHT path (FT, /SEC f ) OH (M/SEC) 

C R (5 ) S FLIGHT PATH ANGLE (RAQ.) 

C S ( 6 ) = HEADING ANGLEfRAO.) 

C S ( 7 ) * ALT I T L'n t (FT,) OR (M) 

C S f 6 ) = X-niSTANC£ (FT.) OR (M) 

C S ( 4 > 3 Y-niRTANCE (FT,) OR (M) 

C S(10)= ACCELERATION ALONG FLIGHT PATH ( F T , / SEC , * *2 ) OR (h/SEC**2) 

C S(U)= TIHE RATE OF CHANGE OF FLIGHT PATH ANGLE (RAD, /SEC , ) 

C 5(12)S TIME RATE OF Change OF HEADING AnG(E (RAD/SEC) 

C 5(13)* RATE OF CLIMB <M # /$FO Cr (H/StC) 

c $(14): SPEED ALONG X-OIRECTT 0 N 

C S(15)= SPFEP ALONG Y-OIRECTION 

C 

1375 S(l) * 6 

1376 $(2) = T ( 2 ) 

i4oo sen = 4,i 

Uflj. S(4) * T(4) 

1403 3(51*0, 

1403 SCSI * 0,0 

1 40 41 S(«> ■ T<5) 

1406 S(91 * 0,0 

1407 CAul INT5(Sr 0 , 2 , l . , t , , 1 ,, 1 , , 1 ,,1 , , DERIVE) 

1420 GO TO 400 

C 

c main airborne integration Loop (2-6) 

c 

1424 ^ NCOL'NT a nCDunT + 1 

1426 ALPHAJ = A( PHA 

1427 THETAJ s THfl AF 

1430 TF{(sC7) * HAPr) ,GF , HMAX)RfcTljRN 

C 

C check PROGRAM PROTECTION LIMITS 

c 

1435 1 F ( AHS ( T ( 5 ) ) ,GT, OSTOPlGn To 9 

144? IF (AHS ( S (9) ) ,GT, D S T D P ) G 0 TO 9 

1446 I F ( S ( 7 ) , l T , - 0 . 1 >Gn T Cl 0 

1450 IF (5(2) , GT , 300, ) bO TO 9 

C 

c obtain atmospheric variables 

C 

1454 300 CAli. AthoS(S{7),OTARS*ANS) 



t/4 

KO 


takoff 

1460 


$ A « A NS ( 4 ) 

1461 


RHQ * ANS ( 3 ) 

1463 


EM a S ( 4 ) /S A 

1464 


HABS * S ( 7 ) - HAPT 


C 

c 

call engine »-tth PWRSET (KENG.O) 

1467 

c 

call Engine (s< 7 ), OTA 88 rEh#PHRSET,RF,KENG) 

1474 


H • » - HfAENP<rS( 3)/3600 t 


C 

c 

bfgtn aerodynamic control 

1500 

c 

K « 1 


c 

c 

CHECK FOR START OF PUL L UP MANEUVER 

1501 

c 

IF ((VEND -SCR) *sqRT(AnS( 7 ) )* YF-LFTR).IT* 

1521 


IF ( ( S ( 7 ) - HAPTT ,LT. HNAN)GO TD 4 

1525 


IF ( IFL V ,FQ. 2 ) GO TO 2 ft 

1527 


IFLY = 2 

1527 


I F ( R T C L ,GT. ROOGO TD 9 

1533 


IFHOCN # EO ■ R ) 


VhArG)GO TO 2T 


1546 


1562 

1565 

1575 

1574 

1575 
1577 
1577 
1601 


1603 

1612 

1620 

1632 

1633 
1 63*1 
1643 

1664 

166*4 

1677 

1700 


C 

C ACCELERATION TO VEnD AT CONSTANT RATE OF CLIMB (28-27) 
C 

*WRITEC6, J040)vEw0 

1040 FORHAT< l*f 29hACCttKRATk TO r.LIHB SPEED OF ,f6,n 
28 IF(5c4)*STNfS(5))*60«0 .IF, RTCL ♦ 1 0 1 > GO TO 29 


C 

C REDUCF alpha TO start acceleration phase 

C 

ALPHA r Al^pHA - 0 A f>T * S ( 3 ) * 0 , 5 

IF(xLF ,LT. 0,85)ALPHA « ALPHA ♦ D ADT*S ( *5 *0 » 25 
00 TO 26 
29 K * 9 

ROC » RTCL 
GO TO 26 

27 IFCJHOIJTE ,CH € 2)G0 TO 41 
XLEMAX c 1.2 

C PULLUP HANFl.'VFO - KTryn PEOUIRFD DADT (27-44) 

C 

CALL PuLLt'P(DADT#KoOFrPWRSET # RENG # HApTiDTABS) 

IFCKOOF *ED, R)GO TO R 
ircTOCN ,FQ* R)WHTTF(6 f 1056)[?A0T 
1056 FORMAT ( ?6H EXECUTE PIH.LUP AT DAD T » if4.2) 

JPOUTE = ? 

GO ™ 4 

41 GO TU<42, 43,44 ) ,*QOt ■ 

42 1F(R(4)*SCRT<A!mS<J))* VELfTR *LT f VEND , AND , S(10) f lT t 

98 

GO TO 4 

43 IF(8(4)ASaRTf aN 8(7) VELfTR .GF, VEND )G0 TO 8 
GO TO 4 

44 IF(SUO) .1 T » 0,02)t;0 TO 6 


, 02 ) GO TO 


C 



TAKn FF 


-P> 

o 


c 

c 

c 

c 

c 

c 

1703 

1706 

1712 

1717 

C 

C 

C 

1717 

1720 

1722 

1725 

17410 
1751 
1753 
1764 
1 7 74 
1775 

2004 

zoos 

2023 

2026 

2036 

2037 

2040 

2041 

2043 

2044 

C 

C 

C 

2044 

2056 

2064 

2065 

C 

C 

c 

2067 

2073 

2074 

2075 
2077 

C 

C 

C 

2100 
210 1 
2111 


m»tASF ANGLE OF ATTACK EACH T I HE BEFORE INTEGRATION STEP TAKEN 
(EXCEPT FOR CONSTANT RATE OF CLl^fi PORTION), FOR LOAD FACTOR 
f LT, 1,0 , additional INCREASE in ALPHA, IF the increase IN alpha 
RESULTS IN any Vim ATTON OF FLIGHT PATH CONSTRA I NTS ( ACCEL. # *LF * THET AF ) , 
ALPHA h'll.l Hf: REDUCED ACCORDINGLY IN DE»IV2, 

4 ALPHA s ALPHA + QAUT*S(3) 

IF(XL* ,LT. 0 , N) ALPHA s ALPHA ♦ DADT*S(3) 

I F ( X L F ,L T , 0 , 6 ) A l PH A * ALPHA F DADT*S(3) 

26 CONTINUE 

roll anG|F control (26-49) 
if(mroli .ec, o ) go to 49 

HEAD = 5(6 ) *57 v 3 
IFOPOLL ,GF, 2)GO TO 45 

OFLPST = ABS(PHI*PhT*Ga0S*CY/(2,*DPHIDT*w*S(4)*C0S(S(5)))) 
JF(AHS(HEaoF-HEAD) ',LT, 57 . 3*Dtl P5l ) GO TU 45 
PHI S PHI +nPHlOTASt‘i) 

TF(AQ5(PHI) ,GT, PHIHAX , AND f PHI f GT , 0 t 0)PHl s PHIBAX 

1F(AHS(PHT) ,GT, PHlHAX , AND , PHI ,LT, 0,0)PHl c *PHIMAX 

GO TO 49 

45 GO 70(46*47*49) * I ROLL 

46 TROLL s 2 

47 1F((ITURN ,r t T # 0 .AND. HEAP .Gt, HE ADF ) ,OR, (ITURN , L T , 0 « AND, 

9 H E A 0 , L T , HE A OF ) ) GO TO 4ft 

PHJ s PHI * OPHJOT#S( 3) 

IFCASS(PHT) p LF., 2./57,3)Ph! s f CHEAOF-HEAO) / A«S ( HE ADF-HEAD)) 

9*2, /57, 3 
GO TO 49 
4H IROLL ~ 3 
PHI * O.C 
5(6) = HLADF/57,3 
HR DLL = 0 
49 CONTINUE 

MAKE InTF.GRATToN STEP 

3 C lNTM(S#6#2,l ll ft,fl,fl t pt,fl,, 0ERIV2) 

I F ( K 99 ) GO tD 9 

T ( 5 ) s S(ft) 

TK ( T 35 .ED, 2)G0 Til 22 

S A VF VALUES FOR nRsTAC|.C HEIGHT INTERPOLATION 

I F ( ( S ( 7 ) - HAPT) » Gf , HDnGo TO 21 

yj * s ( 4 ) 

TJS = 1(51 
MJ e 5(7) 

GO TO 22 

FIND VALUES AT OBSTACLE HEIGHT (21-1015) 


21 135 3 2 

$ 35 = YYYChDt, TJ5,T(5),HJ,S(7)) 

V35 = YYY(HDt, y J , 5 ( 4 ) ^ H J , s f 7 ) ) * $ pR T ( A N 5 ( 7 ) )*0,5920H7 



TAKOfF 

2130 


2151 

2153 


2155 

iitf 

2163 


2200 


22fU 
2205 
2207 
2210 
2211 
221 1 


2213 

2214 
22*2 

2244 

2244 

2250 


2267 

2272 

2276 

2300 

2300 

2302 

2304 


2307 

2307 

2310 

2316 

2340 

2340 

2344 


1FCIDCN , F 0 , 9) 

♦ *RlTfc ( 6rU> 151 S3S, V35 ^ BOCC ^ tlt 

1015 FOR^aT(30^ DISTANCE TO Obstacle height *,F7 # 1#20H SCREEN SPEED (t AS 
91 b,F6 # 1) 

22 IFCIG ,E0. ?)G0 10 25 
1F(IG ,EG, 1 ) GO TO 23 

GEAR REDACTION (22-23) 

IF ( ( S( 7 ) - MAPI) b IT, HGR1G0 TO 25 


|§u s * s s?in ♦ OTGR 
IFUGCK ,E3. 9) 

1 025^ F ORH A T^ 1 27HGF AR^RETWACTION STARTED AT ,F6«l,i7M SEC, COMPLETE AT 
9,F6,l#4k SEC) 

IG * 1 

c gear Drag IhCrE^nT kFOUCFD UnFArLY *lTH time In OTGr StCo*D3 

c 

23 GRCD =CDGFar*(1*0 - ( $ ( 2 ) - TG)/DTGR) 

IFCGPCD ,Gf. O.niGC 10 ?5 
GRCD = Q.c 
IG = 2 
25 CONTINUE 

IF ( JROUTF .CO, 2 1 GO TO 39 
C 

C FLAP RETRACTION (25-16) 

IF (DEI FD >0, C.OJGn TO 16 

10 IFC($(n* ! HAPT) L tlT # HFLAP .OR, S ( 4 ) *$GR T ( AN$ ( 7 ) ) * , 592087 *LT, 

9 VFLAPJGO TO [h 

HfHV E 2 

TIME *(PELEP • »OflFO(IFLAPn/DFL.PDT 
IFtlOfN ,tU. <>) 

♦ ARlTtC6>10iO)'(PFLFP(T' : LAP),TIME 

1050 FORMAT (IX, JRHFLAPS RETRACTEP TO .M.l.RH PEG, IN ,M,1» 5H SEC.) 

15 DELFO » OFLFC - DFlPt>T*$(3) 

IF (DEI FD ,GT. XDELFO(TFLAP))GO TO 16 
DFLED » XDELFPf IFLAP) 

HFLAP = 1 

IftAP a IFLAP ♦ 1 

HFLAP c XHF EAP(IFLAP) 

VFL AP s XVFt AP ( IM. AP) 

c vectored thrust angle deduction ci6-S6) 
c 

16 CONTINUE 

IF(ANGLF ,fcO* 0 # 0)GO to 56 

50 IF ({5(?)-HARn V !tT. HVtCT ,0R, S ( a ) «SORT ( ANS ( 7 ) ) » VELFTR .IT, 

9 WEC 1 ) GO TO 56 
MVECT * 2 

TIME S (ANGLE - XNiJf I VECT )) /L>T VECT 
IFCIOGN .EO. T) 



takoff 


♦WRITE tfa, 10SS) XNilf 7 VECT) * TIME 

toss FORMAT (IX, 3 3H VECTORED THRUST ANGLE REDUCED TO ,F4 t i,9H DEG, IN , 
9F4, 1,5 h SEC.) 

2363 55 ANGLE » ANGLE - DTVECT*S(3) 

IF(ANGLE ,GT, XNucJVfc CT))GO TO 5b 
2372 ANGLE a XNUnVECT) 

237a MVt'CT = 1 

237a IVFC T s IVFCT + \ 

23 7b HVECT a XHVE t T < T Vft T ) 

2490 V VLC T a XVVFC.T ( I VECT) 

C 

C THROTTLE SETTTnG MANAGEMENT (56-39) 

C 

2903 5b CONTINUE 

2403 GO lU(3ir32*33)f JPOwF.R 

24t2 31 IF((S(7) - HAPT) .LT * HPGWFR ,0R f S ( 4 ) * SORT ( ANS ( 7 ) ) * VELf TR »LT , 

9 yPO^ERIGO TO 39 

C 

C DETERMINE PnwFR IwCRtAsF OR DECREASE 
C 

2434 IF(PHRSFT - XPOHERCMPO^ER) )34, 39, 35 

C 

C ADVANCE throttle SETTING LOOP (34-35) 

C 

2441 34 1 POWER a 2 

24<n Tint = 1 On ,M XPO*Fr(mPOwEr) - P W«sE T ) /DjPUP 

2445 SH = XPOWf R (MPONER) *1 00, 

2446 TFCTDCN .FO. 9) 

♦WRITE (fa, li)50)SF T , TIME 

^ 1050 FflRHATc j X t 2BMA0V*NCF THRUTTLF SETTING TO ,Ffa t l,12H PERCENT IN , 

10 9F4.t,SH SEC.) 

2464 32 PWRSF’T a PwRSET + ( D TPUP / 1 0 0 , ) * S ( 3 ) 

2470 If(P*RSFT t l T # XPONf R C HPOWfcR ) )GCI TO 39 

2**73 pWRSFT s Xpn*ER(MPn*fcR) 

2475 MPOWF.R s mPOwT r + 1 

2477 HPDWLP = X M P w R ( M P £ J w E R ) 

2501 VPOWER = XVPHR(MPflwFR) 

2503 I POWER = 1 

2504 GO TO 39 

£ retard throttle setting loop ( 35-39) 

c 

2507 3b I POWER = i 

2507 Tl^E a 100 ,»(P*RsEt - xPOwt R ( MPO w ER )) /DTPDwN 

2513 SET = XPl)WfP(MPnwiFR)*10O, 

2514 IFCIDCN ,EQ* O) 

+V-RITEC6, 1051 )SET, TIME B t 

1051 FORMAT (1X,27HRF. TaPI) THROTTLE SETTING TO ,F5,l f 12H PERCENT IN ,F4,1 
^ f 5H SEC,) 


2532 

33 PWRStT = 

Pv.RSET - (OTPDWN/lOO, )*S(3) 

2536 

If (PWRSFT 

.GT, XPnrfPR(MPOWLR) )G0 TO 

2542 

PWRSET = 

XPn^ER (mpuwf.R) 

2544 

MPOWF R s 

MPO^ER + 1 

2545 

HPOW£R s 

XHPfcR ( Mpn*F«) 

2547 

VPQWf.R z 

X yP^R ( ppflwFR ) 

2552 

1P0WEP a 

1 



TAK off 
2553 


39 CONTINUE 


4* 


C 

c 

c 

2553 

2557 

2576 

2576 

2577 
2600 
2602 


2616 

2620 

2622 

2625 

2627 

2630 

2630 

2637 

2637 

2602 


26*2 

260S 


2607 

2651 

2650 

2660 

2662 

2671 

2672 
2670 


2762 

27cO 

2766 

2766 

2777 
30 IS 

C 

c 

c 

3033 

C 

c 

c 

3030 

3036 

3000 

3000 

3050 

3056 

3^60 

3065 


HEADING CONTROL (39-691 

TRACK a sr}RT<S(*)*S(ft) ♦ S(9)*S(9)) 

IF(($(71- HART! t GE. HHEAO .l)R, TRACK # GE, RANGE)G0 TO 61 
tiO TO 69 
61 MPDU a 1 
TROLL = 1 

HEAOF » XMfc A[) ( iHfc Ain 
IF < I DC K ,FO # 9 ) HR T T F [ 6 * 66 1 HF A OF 
66 FORMAT UX # 22NrtEfi1N TU«N TO HEADING ,F5,1,SH DEG,) 

I wt AO S IHEAD t 1 

RANGE a XPANGE(THFA0)*RNGFTR 

HHfeAO * KHHf ADCIHEAP) 

HEAD a S<M*57,i 
TTURN S 1 

IFtHFAPF t LT. HF:aD)TTURN a -1 
DPhtot s float ctturn)*rpl r at/57 # 3 
69 CONTINUE 

IF(NCDUNT # LT,10)GO to 2 
NtOUNT e 0 

xnADt a (ALPHA - ALPHAJ)/S(3) 

OTHTPT a fTHFTAF - THETAJ)/S{3) 

TA » ESP * THRlfST 
VKTS = S[U)*VELFTP 
FAS a VKTS*SQRT(ANSf 7) ) 

GAMMA a $(51*57,295 
ROC a SC<n*SIN($(5))*hOi 
ROLL = PHI *57.3 
HEAO a S [ 6 } *5 7 • 3 
IFCZDCM ,F0. 9] 

♦ WRITE(6‘, 1002)S(2)#S(6)^S(9),S(7)fVKTS,EAS#EM,S(10) r CL#CD, ALPHA# 
9G A MM A, ROC, XLK, TA , TmFTAF, ROLL , HEAD 

I PAGE = IPAGE *1 

ifcipagf ,lt # ? -page. > go to ^ 

IPAGC a 0 
IFdOC.N ,F:K 9) 

♦ MRITF<6,99tn 

IF(IUCN ,FP # 9 .AND, tUNir .E&, 1) HR 1 TE ( 6 , 1 00 0 ) 

T F ( 1 0 C N m ru\ 9 .AN.), I UNI T .FO. 2 ) HR I TE ( 6 , 200 0 ) 

END OF MAIN AIPRORnF INTEGRATION (OOP, 

GO TO '2 

CONVERSIONS AND PRINT (JUT. 

ft TA a F nR * THRUST 
yKyS * SCO)*vELFtr 
EAS a VKT5*S(JRT (AMS(7) ) 

IFlKODf ,Nf, 3 ) F A S a VEND 

GAMMA a $(51*57,295 

ROC « Sf«)*SlM(S(S))*60 # 

ROLL a PH 1*57 ,3 
HEAD a S(6)*57 t 3 



TAKOFf 


3067 

IF(IDCn 

t 


♦WRITE ( 6, 

! 


9GAMMA f ROC 

3155 

WRITfc (*, 

1 

3161 

If CKQOE 


3176 

If (KODE 

f 


F 0 « *) 

002)SC?),S(tn,$(9),sm,VKTS,EAS,E:M,8(i0> f CL # CD # ALPHAf 
,XlF # TA,THETAF.ROLl ,HEAD 


052) 


J 'I uiU T TP t A . 1 




1052 FORMAT*/, p(, J4HFN0 OF TAKEOFF) 

1053 FORMATnX^6HTH«aTTt iMti REQUIRED TO MAINTIAN CONSTANT SPEED) 

1054 F0rMAT(1X,59HDE5;IrED END S^ED NOT ATTAINABLE At SPECIFIED POWER 
RETTING) 

3210 RETURN 

321 1 9 CONTINUE 

3211 IF (R TCL ,GT. ROC ) WR I TF C 6 , 996 ) 

996 FORMAT (/,2X,66HCAnN 0T ACCEL , AT INPUT H/C fRTCL), TRY VALUE ,LT t 
9 AST R/C PRINTED) 

3223 WRITE (6 , 997 ) 

991 FORMAT (/, 1 X,36H*** ABNORMAL TERMINATION OF TAKOFF ***) 

3233 RETURN 

3233 FNt> 


S 


L 


SUBPROGRAM length 
04566 




FUNCTION ASSIGNMENTS 
YYY - 000020 


STATEMENT ASSIGNMENTS 


1 


000770 

2 

m 

001425 

7 


0C0150 

6 

m 

003f»35 

21 


002101 

22 

m 

002152 

26 


001547 

29 

* 

001575 

35 


002506 

39 

m 

002554 

45 


001776 

46 

m 

002005 

55 


002364 

56 

w* 

002404 

102 


001057 

103 

m 

001050 

201 


003454 

202 

m 

003466 

207 


003600 

206 

m 

00360 7 

213 


0Q3644 

214 

m 

003656 

300 


001455 

5 00 


003434 

999 


003677 

1000 


003753 

1025 


004130 

1030 

* 

004143 

1053 


004225 

1054 

w 

004234 

BLOCK 

NAMES AND LENGTHS 



UNI V 

* 

0O0030/01 

AFRO 

- 

000020/02 

UNIT 

- 

000002/07 





3 

p 

002045 

4 

- 001704 

9 

p 

003212 

10 

- 002223 

23 

m 

00220? 

26 

- 002212 

31 

m 

002413 

32 

- 002465 

41 

m 

001635 

42 

p. 001644 

47 

p 

002006 

ad 

- 002040 

61 

m 

002577 

66 

• 004210 

106 

"i 

001000 

no 

- 001061 

203 

m 

003504 

204 

• 0 035 1 3 

209 

m 

U 036 16 

210 

* 003622 

215 

m 

003663 

216 

• 003667 

501 

m 

000102 

996 

- 004244 

1002 

p 

004045 

1009 

- 003710 

1 040 

P 

004066 

1050 

p 004164 

1055 

m 

004152 

1056 

» 004102 

ROLL 

• 

000002/03 

f XCMNG 

p 000004/04 


5 

* 000037 

6 

- 000127 

15 

* 002270 

16 

• 002310 

26 

« 001720 

27 

• 001*00 

33 

♦ 002533 

34 

- 002440 

43 

* 001665 

44 

• 001701 

49 

• 002045 

50 

- 002317 

69 

» 002640 

101 

- 001013 

113 

• 001134 

120 

• 001312 

205 

. 003531 

206 

- 003557 

211 

* 003626 

212 

- 003*35 

217 

- 003673 

299 

• 0003*0 

997 

• 004255 

996 

> 003733 

1010 

- 003740 

1015 

- 004117 

1051 

p 004176 

1052 

■ 00422J 

1057 

p 003725 

2000 

• 004010 

XHOLL 

p 000041/05 

XFLATE 

• 000121/0* 


VARIABLE 

ALPHA p 

ASSIGNMENTS 

mm/ oi 

i 

ALPHA J 

- 00452? 

ANGLE 

m 

000005/02 

ANS 


004323 

CDGFTP « 

CL 

p 000013/02 

CY 

w 

000012/02 

DADT 


004454 

DFLPD7 - 

0O443G 

DPHTDT 

- 004547 

DS TOP 

p 

004461 

DTA0S 


004431 

DTPDWn - 

004433 

DTPUP 

• 004434 

DTVECT 

m 

004435 

EAS 


004525 

ENP 

000021/01 

enpout 

p 00453? 

E YEW 


000020/01 

E YEwNU 


004436 


CD 

DELED 

DTGR 

EASTO 

G 


000014/02 

000006/02 

004432 

004531 

000001/07 


COGEAR 

DELPSI 

DTMTDT 

EM 

GAMMA 


004427 

004946 

004927 

000010/01 

004562 



TAX OFF 
GRCD 

000016/02 

HASS 


000002/02 

HE LAP • 

004534 

H&R 


0044U0 

HPOEIER - 

C04536 

HVECT 


004540 

IEL Y 

004473 

TG 


004504 

IROLI - 

604505 

I TURN 


004551 

JJl 

004467 

J J2 


00^470 

METRIC - 

004453 

mFLAP 


004476 

NCOUNT - 

004521 

NEO 


000036/05 

PWRSET - 

0 0 45 1 S 

os 


000001/02 

ROC Ml N - 

000001/03 

roll 


004563 

SA 

004520 

SET 


00456Q 

T A 

000022/01 

TEMP 


00451 ? 

THETAS - 

004523 

thrust 


OOO 003/02 

TJ5 

004563 

tow 


004514 

veletr - 

004457 

VEND 


0 0 0 0 L Cl 

VMARG • 

004533 

vpnwEH 


004537 

WG 

000023/01 

was 


004513 

XHHEAD 

004415 

XHPWP 


004352 

XPOwER - 

004164 

XRANGE 


00441 0 

ZERO - 004516 

START OF CPNS1AMS- 
7600 COMP It A T I OK "• 

003423 
RUN 76 

TMHS- 
Lt'VEt 9H 


ROUTINE COMPILES IN 060200 


MART 


004437 

HOT 


004452 

HHEAD 


004543 

HJ 


004554 

I 


004464 

lOOhN 


004472 

IHEAD 


004507 

TOUT 


004442 

tuntt 


000000/07 

TUP 


004471 

JHOUTt 


0044 76 

K 


000120/06 

wpDWtrt 


004500 

MROLL 


004506 

NT A&E 


004456 

PHI 


000000/03 

range 


004542 

RHO 


000015/02 

ROLL MX 


004450 

ROLHAT 


004451 

SWOLl 


UOOOOO/04 

SH 


000015/01 

TG 


004555 

TGU 


004556 

THTF1 Y 


004446 

THTSCP 


004447 

TRACK 


004561 

UM 


004443 

Vpl AP 


004535 

VJ 


004552 

VVECT 


004541 

V35 


000002/04 

XDADT 


004526 

XOtLFD 


004333 

XMVFC T 


004376 

*LE 


0001 16/06 

XVFL AP 

m 

004145 

XVPWR 


004357 


00 42 1 5 I wo I RFC TS« 0 04321 

74/07/15, 


HEAD 


004545 

MEADE 


004550 

HH AN 


004441 

HMAX 


004455 

IE AST 


000017/02 

IEL AP 


004474 

IPAGE 


004465 

I POWER 


004477 

IVECT 


004502 

135 


004503 

KENG 


004466 

KODf 


004544 

MU 


000037/06 

HVECT 


004501 

PH I M A X 


004610 

PHARG 


004444 

rngftr 


004460 

ROC 


0001 13/06 

RTCL 


004445 

S 


000000/06 

535 


000001/04 

T 


000000/05 

THEMAX 


000114/06 

THETAF 


000115/06 

TIME 


004557 

TINT 


004511 

UNITS 


004463 

VEL 


000000/02 

VXTO 


004530 

VKTS 


004524 

w 


000006/01 

WF 


000007/01 

XHEAP 


004422 

xhflap 


004340 

XLFHAX 


000117/06 

XNU 


004371 

XVVECT 


004403 

z 


004517 



DEPJV1 

SU0ROUTINF OtPIVi 

c 

C SUBROUTINE OFRIVi COMPUTES THF ACCELERATION Tf6) Fqr THE GROUND ROLL 

c 

PEAL HU 

COMMON /IJNIV/ NPC ,N$C f IOC 9 H f ST rR r 

UF ,fH , VMO ,£HMO f ALPHLO,CLAtPH,SW , AH ,8 , 

?ETtW # F NP .Ta ,J*GS * KWR I Tfc , 01MC 4 

3,HSIZE 

COMMON /*£*"/ VEL*QSfHAB5fTHHUST,TVECT^ ANGl fc * DELF D f OELSPL , ALPHA , 
9rV,CY>Cl.#tp,RH0,(;PCD, IFA5T 
COMMON /XPOLL / T (30) # NEO# MU r nW£ V 
COMMON /UNIT/IilNlT^O 
OS c 0 , 5*PhnASrt*T f 4 ) *T ( 4 ) 

0 t r COS ,FQ . O.lUS = 0.1 

7 IFAST S 1 

10 CALL AHDpy K 

11 T(6) * (G/kO*(-k*MU + US*<CY*MU - CXH 

iiO T(7) = TC«) 

rM RFTUHN 

*2 END 


SUBPROGRAM LENGTH 
00040 

FUNCTION ASSIGNMENTS 
STATEMENT ASSIGNMENTS 
BLOCK NaM£S AND ItNGlHS 


UNIV 

000030/ A 1 A£PO 

- 0000P0/02 

xroll 

- 000041/03 

UNIT 

• Q 0 0 00 2/ 04 

variable 

assignments 






cx 

ooooii/o? cy 

• onooiz/na 

G 

- 000001/04 

IFAST 

* 0000 1 7/0^ MU • 000037/03 QS * 000001/02 

RHO 

00001S/02 SW 

- oooois/ot 

T 

- 000000/03 

W 

- 000006/01 

START OF 

CnNSTAMS<-000025 

TEMPS-- 

OOOOJO 

INOTRECTS- 

000040 



7600 COMPILATION — RUN76 Lf VH QH 74/07/t5 # 

ROUTINE. COMPILES in 0 4400 0 



DER1V2 


SUBROUTINE OfcftlVZ 

C SUBROUTINE n*RrV£ COMPUTES the TIME DERIVATIVES FOR THE AIRBORNE 
C PORTION OF THE TAKFOFF AND MANAGES THE FLIGHT PATH CONTROL, 

C 

COMMON /untv/ NPC .NSC 1 1 DC ,H ,$T t R ♦ * » 

twF i EM . V MO ,EMM0 , ALPHLO,CLALPHrSW ,AR ,B i 

PF.YEW f ENP . T A ,w G ,*GS , KNRI TE . DLMC4 

3.KSIZE 

COMMON /aero/ VELrQS'HABS'THRIJST'TVECTr A NGLE * PELF D r DELSPL # ALPHA , 
9CX.CY.Cl, cn, RHn,GRCn, IF A ST 

COMMON /xF LATE/ SO l 5)rRnCrTHEMAK,THETAF.X!>,XLFMAX#K 

COMMON /ROLL/ PHTrPOCMIN 
COMMON /UNJ T/I!IN(T, G 
MER = 1 

2 QS = O.S*RHn*S^*S(4)*S( / 4) 

5 S'KTS - S(a>*0*592097 

7 tFCTUMT # F«, 2)VKTS = R C«) 

C 

C CONSTANT RATE tjf r l I MM PORTION 

c 

15 1 IF ( K ,Fq, 9ICAII CLtMB(R0C.S(5)f VKTSf NER3 

20 1 F ( NE R , NF , nHRITE(6 J r6b6>ALPHA,R0C 

66b FORHATn x ,30H*AFRROR IN CLIMB - ALPHA , ROC »,2F10,2) 

52 26 I F A 5T = 0 

55 CALL ARQDYN 

C 

C CHECK FUSELAGE AnGIF, IF THF.TAF ,GT, THFMax, REDUCE ALPHA (26W62) 
C 

3« 261 IHLTAF * s(S)*57,?9S + Al,PHA - E Y E W 

a 0 IFCTHETAE ,1F, THFMAXIGO TO 262 

<♦2 ALPHA * THEMAX ♦ E/Ffc - SC5)*5?,295 

a5 GO TO 26 

C 

C CHECK load FACTOR, IF xLF # GT, xLFMAX. REDUCE ALpHA (262-263) 

c 

as 262 XI. F 5 COS*CV)/^ 

a 7 IF(XLF ,LE. xLFmax)Go to 265 

S2 ALPHA 5 A | PHA - 0.05 

5<l GO TO 26 

c 

C CHECK ACLFLFkAtToN, IF Stm.tT, 0,0, REDUCE ALPHA (265*50) 

C IF TN CONSTANT RATE CF CLIMB, RFDUCF ROLL ANGLE (MIN OF 5,0) 

C 


55 

263 

s C 10 ) = (G/W)* (~CX*OS “ W*SlrjfS(5) )) 

6a 


IF($(1 0) ,GF , 0 , O ) GO TO 30 

66 


I F ( K ,EO, 9)G0 TO 26« 

IS 

265 

WS CpB* Al .rv: TO 99 

7a 


GO TO 26 

75 

26a 

IF(ABStPHT) .1 T # S.0/57, 3)Gt) TO 265 

l o i 


IF (PHl ,GT, U,)PHT * PUT - 0,1/57,3 

1 03 


IF(PHI ,LT 4 0,)PHT r PHI t 0,1/57,3 

t 06 


Gn TO 1 

107 

30 

CONTINUE 

1 07 


AWG a w*CPS<$(5))/rCY*QS) 

115 


I F ( ARG , G T , 1,00) ARG * 1,00 



-t* 

00 


0ERIV2 

121 

142 

154 


164 

167 

173 

174 
206 
211 
217 
226 
226 

235 

236 


IF(fttS(PHT) # GT , 0 # ,AND # S (4 ) *SlN( S (5) ) *60 t ,lT, ROCMlNjPHl ■ 
+ CPHl/AB5CPHin*AC0Sf ARC 5 

SC11) * (G/(**Sf4) ) )*fCY*0S*C0SfPHn ** W»C0&CS(5))) 
IF(S(in*57,3 ,GE. *1.0360 70 40 


IPCPHI ,tq. 0.)Gn TO 40 #c , , 

IFCPHI # LT , 0 f ) Phi a PHI + 0 # 1/57 1 3 
IF (PHI ,GT . Ot)PHT * PHI - 0 f 1/57 « 3 
IFUBSCPwn .IT. 0. 15/57,3)PHI * 0 t 
GO TO 30 

40 $ ( 12) s (G/ew*S(4)*C0StS(5))l ) * ( C Y*QS * S I N ( PH I ) ) 

S ( 13) s S(4)*SIN{$(5) ) 

$ f 1 4) b S(4)*C0S(S(5))*C0S(S(6)) 

SU5) * 3fw)*C03(S(5))*SlKi(S(6)) 

RETURN 

99 „R]TE(6,66)S(&) 

66 formaT(1x,46h*** unable to maintain accel. ,ge # o,o# 

K * 45 
END 


DV/OT*f F9, 5) 


SUBPROGRAM LENGTH 
00347 

FUNCTION ASSIGNMENTS 


statement ASSIGNMENTS 

1 - 000014 26 

261 - 00003b 262 


000033 

000046 


2^3 


OOQf 10 
0000*6 


% 


* 000175 

* 000076 


Sis 


mm xu 


BLOCK NAMES AND LENGTHS 

UNI V - 00003n/ni AERO - 000020/02 XFLATE • 000121/03 ROLL 


000002/04 UNIT 


000002/05 


VARIABLE 
alpha - 

ASSIGNMENTS 

00001C/02 

* 0 n 0 346 

ex 

000011/02 

CY 

IFAST * 

000017/02 I U N I T 

- 000000/05 

K 

000120/01 

NER 

PHC - 

000015/02 «QC 

- 000113/03 

ROCMIN * 

000001/04 

S 

thetaf - 

0O0IJ5/O3 VKTS 

- 000345 

W 

OOOOOG/OI 

XLF 


000012/02 

000344 

000000/03 

000116/03 


EYE W • 000020/01 G 

PHI * 000000/04 QS 

S h - 000015/01 THEMAX 

XLFMAX • 000117/03 


START nr CON ST ANTS* 00024) 


Tt Hp$..*000$0O 1 NDI NEC TS^OO 0 336 


7600 CriMPll AT I ON — *IJn 76 LEVEL 9H 74/n7/15 # 


MiiS 


000001/05 

000001/02 

0001H/03 


ROUTINE COHPItfS IN 004300 



CLIHB 


SUBROUTINF CLlMRfROC pGAHHA, VKT5 # NER) 




c 

c 

c 

c 

c 

c 

c 

c 


6 

7 

11 

n 

12 
14 
17 
23 
25 

30 

31 

32 
34 
*0 
ft 6 

« 

104 

104 

105 
loft 
107 


subroutine ci finds the required alpha to Ely at the constant 
rate Of CLIHp RTCL, GIVEN THE THRUST AND VELOCITY. CHANGES IN FLIGHT 
PATH ANGLE ARE FaJRLY INSENSITIVE TO VARIATIONS IN ANGLE OE ATTACK, 

AS A RESULT, THE COMPUTED RATE Of CLIMB WILL DIFEER SOME {USUALLY LOWER) 
THAN THE DESIRED VALUE rtcl, 
subroutine zerjvb is a ZERO»F I NDER , 

RE A l NU 

COMMON /UhIV/ NPC 9 MSC 9 IDC fH r ST ,R f W , 

1 wF , EM # VMQ f f N M 0 , ALPHLO,CLALRH r sw , AR #8 9 

HfcYEW ,ENP ,ta »WG #WGS ,KWRITE,DLMC4 

3#KS IZE 

COMMON /aero/ VELr05#HABS,THRUST,TVECT/ ANGLE , DELE D t DEL9PL * ALPHA , 
9CXfCY#CLfCD*RH0,GRCD/ TF AST 
COMMON /ROLt / PHI, ROC HI N 

COMMON /umit/iunit,g 

NER * 1 

TOL = 0,01 

STEP c 1.0 

TX e 0 

JC » 0 

FACTOR r 2.B561 

IFdUNIT *, FO* 2 ) f AC T OR * J.Oo 

QS * 0,5*RHn*VKTS*VKTS*SW»FACT0R 

raunu z PQQ • 

40 IFCJX .EG', o ) F'RRM i s ERROR 
1 F AST * 0 
CALI arodyn 
ALPHX ■ ALPHA*, 0174*55 

ERROR a OS* ( C YaCO* { GAMM A ) « C OS ( PH I ) " CX* S I N ( G AMM A ) ) • W 
IF ( ABS(ERROR) .LT, 0,0025)GP TO ftO 
I F ( JX , FQ. 2 ) GO TO 60 

^F , ?ltC ZE Gf yR ^?GS R fo R 65 1 # ALf>HA ' STE ' P ' T ° L ' JCr JX) 

GO TO 40 
60 RETURN 
65 NER = 9 
RETURN ■ 

END 


SUBPROGRAM LENGTH 
00147 

function assignments 

STATEMENT ASSIGNMENTS 


40 

- 00002* 60 

* 0001 OH 

65 

- 000106 


BLOCK 

UNIV 

NAMES AnO LENGTHS 
- 000030/01 AEP(1 

- 000020/02 

RUU 

• 000002/03 UNIT 

* 000002/04 



cu*e 


. 000012/02 ERRB1 • 000145 ERROR 

• 0 0 0 1 *• 2 JX a OOOI 41 NU 

a 000140 SW a 000019/01 tOL 

START OF CPNSTAfJTS-OOOl la TEMp$-*0 0 0 1 ?2 I MO I RECTS*000 1 Jfe 

7600 COMPlt.ATlOix -- RUN76 I.EVtl RH 74/07/15, 

ROUTINE COnPILfS IN 044000 


VARIABLE 

ASSIGNMENTS 



• 000011/02 

CY 

ALPHA - 

00001(5/02 

ALPH* 

- 000146 

cx 

FACTOR - 

ooom 

I FAST 

- 000017/02 

TUN1T 

• oooooo/oa 

JC 

PHI * 

w 

OOOOOO/fli 

oooooe/oi 

os 

- 000001/02 

RHO 

• 000015/02 

STEP 


• 000144 

• 000114 

• 0001)7 



PULLUP 

SUB ROUT INF PUl LUP ( DADT, KODE# PWRSET , KENG# VEND , N APT, DTABS) 

C 

C SUQROUTlNF PULLUP DTERMInE$ THE TIME RATE OF CHANGE OF THE ANGLE OF 
C ATTACK DAOT REQUIRED TO BRING THE AIRCRAFT FROM THE CONSTANT RATE OF 
C C t- 1 M 9 TO ThF FINAL CLtMB SPEfD(t,F, REDUCE ACCELERATION ALONG FLIGHT 

c path to zero and thus increasing flight path angle:, the pullup 

e MANEUVER STARTS AT SPEtD OF VEND - PHARG*VEND<E,G, VEND«250 t # PHARG*0 1 04# 
c • START OF PULLUP AT 240, KNOTS ) , THE SUBROUTINE DOES THE VlRY SAME 
C INTEGRATION LOOP as TAKOFF, WITHOUT ANY print OUT t AND VARIES THE 

C VALUE OF DAOT UNT It PROPER VALUE FOUND, THE SEARCH FOR THE REQUIRED 

C DAOT IS FOUND BY A BISECTION TECHNIQUE, 

C 

t IF puLLUP FAILS HlTH DA0 Tb4 § 0 (MAY, VALUE ALLOWED)# USER SHOULO INPUT 
C LARGE R VALUE for PHARR, THF PROGRAM will generally over-shoot the 
C END speed by A KNOT or so, 
c 

EXTERNAL 0ERIV2 

COMMON /t*lv/ NPC #NSC , I DC ,H #ST ,R , W t 

lwf #Fm #VMO t EMMO # ALF'HLO^CLALPHfSW # AR ,B # 

?EYEW #FNP ,TA #«G # WGS $ KhRIT E# DLMC4 

3# KSIZE 

COMMON /AFRO/ vel#gs#habs#thrust# TVFCT, ANGLE#DELFD*DEL9PL#ALPHA, 
9CX#CY#CL#CD#RH0#C.RCD# I FAST 
common /XFLATE/ Sf ;S),ROC#THFMAX#THFTAF#KLF#XlFMAX#K 
common /Un IT/I UN IT, G 
0 1 mF NS I ON tKSiH) 

VELFTfi r 0,59?&B7 


12 


IFUUNIT .EG, 2)VELFTR s 1,00 

16 


S2J * S ( 2 ) 

17 


S4J = S(4) 

21 


S5J 5 S(5I 

22 


S6J s S f 6 ) 

29 


S7J = S C 7) 

25 


S8J * SC8) 

27 


S RJ * SC9) 

30 


WJ * w 

32 


ALPHA J 5 ALPHA 

33 


kodf ■ 1 • 

34 


ILOOP = o 

35 


DADT = 4,0 

36 


OADTIC = U.O 

37 


ITIME a 1 

40 

50 

ILOOP * II OOP T 1 

42 


IF ( ILOOP ,GT. 15) GO T 0 95 

45 


ISTAPT r ) 

46 

5 1 

SC2) • S2J 

47 


S ( 4 ) r S4J 

51 


S(5) = S5J 

52 


3(6) » 96J 

54 


S (7 ) = S7J 

55 


see) * sb j 

57 


s(9) » 391 

60 


m a w.J 

62 


ALPHA s alphaj 

64 

300 

CALL ATM0S(S(7),DTAH5,ANS) 

75 


3 A a ANS ( 4 ) 



cn 

to 


PULUJP 

76 

77 
101 
103 

112 

131 

140 

m 

ms 

146 

ISO 

ISO 

161 

173 


200 

207 

220 

220 

222 

215 

236 

241 

252 

255 

255 

256 

256 

257 
262 


262 

263 

264 

265 

266 
267 

277 

300 

300 

3 1 0 
31! 
3U 


RHO * ANS(3) 

EM a S ( 4 ) / 5 A 
HABS * S ( 7) - HAPT 

CALI ENGINE (St 7),0TA63, EM, PWR5E T , *F , KENG ) 

XE(ISTART 9 ME 4 21 

OC ALL lNTS(S.6»2.l t #U#l«rl.f r)ERlV2) 

IF ( 1ST ART .EQ , 3) RETURN 
I ST 3 2 

* * w - wf*ENP*St3)/3600 f 
K • t 

ALPHA = ALPHA ♦ OAOT*S(3) 

IFfXLF ,LT. 0*9)ALPHA s ALPHA t DADT*S(3) 

IF(XLF ,IT, 0,G)ALPHA a ALPHA ♦ DADT*S(3) 

CALL MTmcS,6,2, 1 m I 1 « *1 »#1 " 1 ,* DER I V2 ) 

EAS r S(4) aSQRTC ANS(/n*VUFTR 

c TEST Fuft VARIOUS MO CONDITIONS ( -110) 

c 

IF(S(U).n, 0.02 ) Gfi TO 100 

IF(FAS t GF. VEND t 0.5 .AND. IT I HE .FQ, l)GO TO 100 
GO TO 300 

MO IF ( IT I HE .KG. 2)GD TO 101 

IF(S(10)*GT. 0,02 .AND. EAS .GE, VEND ♦ 0.5)GO TO 150 

I T I M E « 2 

lOl IF(DaOT .IT. n. 03)GO TO 160 

IF (EAS f GT. VEND .AND, FAS . LT t VEND + l t O)GO TO 200 
IF(EAS .lt. VEND1G0 TO 110 
IFfDADT .ED. 4.0) GO TO 99 
DADTLO 3 DAM 
GO TO US 
HO DADTUP * OAPT 
115 OADT a 0.5*(PADTUP ¥ DADTLO) 

GP TO 50 

C *qOE a 2 ENn spKEP REACHED, fllJJ TOO *UCM THRUST AVAILABLE - THROTTLING 

c KO0E * F 5 R -°C*^OT peach DESIRED E*D SPEED at SPEtIFIED PU*ER SETTING - 

c PI.'ULUP DONE TP. ZERO ACCELERATION 
c KOOE a 9 - PROGRAM FAILS 
C 

ISO KODP - 2 
GO TO 200 
160 KODE * 3 
200 1ST ART t 3 
GO TO 51 

99 ‘aiPITL ( 6 # 9 A ) , _ Ar , , 

9fc> FORMAT (tX,4 y H*** FAILED IN PUILUP * TRT INPUT PhARq a 0,1) 

KPDE a 9 
RETURN 

95 W»ITE<6,96) 

96 F0RHaT( 1X,?9H*** FXCESStVE LOOPS IN PULLUP) 

K OOF * 9 

RETURN 

END 



RULLUP 


subprogram length 
00406 

FUNCTION ASSIGNMENTS 
statement ASSIGNMENTS 


50 

- 000041 

51 

- 

00004? 

95 

m 

000101 

96 

* 

000340 

96 

. 

000331 

99 

* 

000270 

100 

- 000221 

101 

- 

000237 

110 

•t 

000257 

us 

V 

000260 

150 

m 

000263 

160 

«r 

00026S 

200 

- 000266 

300 

- 

000065 













BLOCK 

NAMES AnP LENGTHS 















UNI V 

- 000030/01 

AERO 

- 

000020/02 

XHATE 

m 

000121/03 

unit 

w 

000002/04 







VARIABLE ASSIGNMENTS 















ALPHA 

• 000010/02 

alphaj 


000376 

ANS 

m 

000354 

DADTLQ 

m 

000400 

DADTUP 

m 

000405 

OTA05 

• 

000000 

FA8 

* 000404 

Em 

* 

000010/01 

l NP 

m 

000021/01 

FACTOR 

- 

000365 

HABS 

m 

000002/02 

ILOOP 

• 

000377 

istart 

* 00040? 

ITT ME 

m 

000401 

IUNIT 

m 

000000/04 

K 

m 

000120/03 

RHO 

• 

.000015/02 

S 

m 

000000/03 

SA 

- 000403 

52 J 

m 

C 0 0366 

$4J 

m 

000367 

S5J 

P 

000370 

S6J 

m 

000371 

S7J 

m 

000372 

36J 

* 000373 

S9J 

m 

000374 

VELETR 

m 

000364 

w 

* 

000006/01 

NF 

m 

000007/01 

HJ 

w 

000375 


KlF - 000116/03 

START OF CONST ANTS-000310 T t MPS--00 Q 34 5 INDIRECT 3*000353 

7600 COMP 1 1 AT I ON *- R|JN7* LtVH 9B 74/07/15. 

ROUTINE COMPILES IN 044600 
Ui 



«RJVB 


SUBROUTINE ZEpJVB (ERROR# ERRH t , OR I VfcR r STEP t TOL , JC # JX ) 
ERR * ERROR 


Ln 


il 


IF(JC f GT t oinn TO 10 


14 


JM»0 


14 


JP*0 


15 


JF»0 


15 


JX«0 


16 


JABO 


17 

to 

JCB.TC + 1 


20 


IF ( JP,GT 1 0 ) GO TO 20 


23 

12 

JP*JP*1 


25 


ORM1BOHIVFR 


25 


DRIVER » DRIVFR ♦ STEP 


27 

15 

RETURN 


30 

20 

CONTINUE 


30 


IF (.IF *GT . 0 ) GO TO 45 


33 


IFCERRM1 , 1 T # 0. .ANO f FKR.GTiO, )GO 

41 


IF(FRRHl g GT # 0. f ANO # ERR # LT t O,)GO 

47 


IF (ERR, L T . O , ) 60 TO 25 


50 


IF(JM,GT # 0)GO TO 22 


52 


TF(FRR # GT,FRRHUGO TO 22 


54 


GO TO 12 


W 

22 



56 


DRIVER 3 BU - STEP 


57 


JMBJM+I 


61 


GO TO IS 


61 

25 

TF(Jh f CT,0)G0 TO 22 


64 


lF(ERR f LT*TRP w nGP TO 22 


66 


GO *0 12 


66 

30 

IFUJRN1 f GT.r>Rt VER1G0 TP 35 


7? 


6L«ORHi 


7? 


HIJbORT VfR 


74 


GO TO 40 


74 

35 

RUbORMI 


75 


bLBpRT VER 


76 


FRRmJ x ERROR 


100 

40 

JX B 1 


1C l 


I F ( Jf ,GT. 0 ) GO TO 45 


104 


JF * JF 4 1 


105 


DRIVER a Bl * 0.5 * (Hu - Hi ) 


no 


RETURN 


in 

45 

IF(FRPOR * EPRMi ,LF* 0,0) GO 

10 

113 


BL ■ DRTVFR 


114 


ERRM1 s ERROR 


115 


GO TO 47 


116 

46 

HU * DRIvER 


117 

47 

DRIVER s PL ♦ 4.5 * ( HU -BL ) 


122 


IFCABSCQU - bl.) # n t THDJX = 

2 

127 


RETURN 


130 


END 



SUBPROGRAM LENGTH 



ZEPJVB 

00153 

FUNCTION ASSIGNMENTS 
STATEMENT ASSIGNMENTS 


10 - 000020 12 

• 0 0 0 0 ? 4 

15 

• 

000030 

?0 

- 000031 

22 

• ooooss 

25 

• 000062 

SO > 000067 3b 

• 000075 

40 

m 

000101 

45 

• 000112 

46 

* 0001 17 

47 

• 000120 

BLOCK MAKES «ND LENGTHS 











VARIABLE ASSIGNMENTS 
BL - 000152 UU 

- 00M5J 

ORM| 

• 

000)50 

ERR 

• 000143 

JA 

• 000147 

Jf 

• 000146 

JM - 000140 JP 

- 000145 

JX 

m 

000000 








START OF C0NSTANTS-G00131 Tt MPS—OO 0 1 35 I NO I REC T5*000 1 A3 

TbOO COMPUATinN — RUN76 LEVEL 9b 74/07/15, 


ROUTINE COMPILES IN 044000 



ATMOS 


t/i 

On 


C 

C 

C 

c 

C 

t 

t 

c 

c 

c 

5 

51 
la 
15 
21 
26 

33 

43 

45 

46 
SO 

52 

ll 

74 

75 
77 

100 

103 

105 

107 

110 


SUBROUTINE ATMUS (HH f OTABS, ANSI 
COMMON /U^IT/IUNI T,G 
DIMENSION AMStfl) 

WH-A( TITUOe IN FEET 

0UB5*TEMPFRMIW TNCRFmEnT FROM STANDARD TEMPERATURE 
ANSCn-TEMPERATURE (RANKTNE) 

AN$(?)~ PRESS URL (PSF) 

AnS(3)«0ENS!TY ( 5LUG/F1 3 ) 

ANS(4)-SPEF0 CF SOUND Tn FT f /SEC # 

ANS(5)-KINEMATIC VISCOSITY (FT2/9EC) 

ANS(6)-PRfcSS0RE RATIO 
AN9(7)-DENSITY ratio 
AN$( fl)«.TEMPEPATURfc RATIO 
HJ b HH 

IMIuNIt .fG. 2)HJ S HJ / 0,3040 
CTABJ c OTABS 

IFUUNIT ,EP. 2 ) D T A B >? S 0TABJ*(9 f /5,) 

TNETA a l ,->OoOO607‘i*MJ * DTA6J/S1B,67 

OFtTA s ( l # - # OOOOn6A75*RJ)**5 f 2561 

IF CHj.LE, 360*9,) GO TO 4 

THETA = .7519 t t)TABJ/ 5 ia ,67 

DELTA * ,2?336*f XP(C36OS9 t -Hj>/2O706,) 

4 STGMA = DELTA/THtTA 

ANStn = THFTA ft 510,67 

ANSC 2 ) * DELTA ft ?J 16,22 

A w s ( 3 ) s SIGMA * .0023769 

AN 5 ( 4 ) * 1117.061 *SDRT (THETA ) 

P * AN$f2WW]U 

ANS(S) = , 27 o 550 E-O 6 *ANSCl)*SORT(ANSCU)/ , (PAO,tl 96 ( 72 /ANS{l))) 

A N $ f 6 ) * OF I TA 

ANSf7) 5 SIGMA 

ANSC0) * THETA 

1 F ( I jjN| T t F_ Q . 1 ) RETURN 

ANSC3) * ANS{3)ft5tS,3S 

A N S C 4 ) = ANS( 4 ) * 0,3040 

RETURN 

END 


subprogram length 
00172 

FUNCTION ASSIGNMENTS 

STATEMENT ASSIGNMENTS 
4 * 000044 

Block names and lengths 
UNIT - 000002/M 

HAP14&LE ASSIGNMENTS 

DELTA - 000167 DTABJ - 000165 HJ « 0O0164 IUNIT . 000000/01 P - 000171 SIGMA - 000170 

THETA - 000166 



